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ABSTRACT 


i 


The theory of the computable a-posteriori error estimate for a finite element 
method is developed. Among other things, it is shown that the error estimate is 
very reliable and the ratio (called effectivity index) between the estimator and 
the true error approaches one. Numerical examples computed by program FEARS 
(Finite Element Adaptive Research Solver) of the University of Maryland, illustrate 
the effectivity and reliability of the estimators. 


A 





1. INTRODUCTION 


Recon! Iv an i nr runs i ng into rent in the finite element computations is 
being focused on the reliability of the results and the quality of the used 
meshes and elements. 

During recent years at the University of Maryland, the studies were under¬ 
taken which focused toward the development of a finite element system having 
the following features. 

a) The solver supplies the user with a reliable and accurate information 
about achieved accuracy in the desired norm. 

b) The solver constructs adaptively meshes which are leading to the highest 
possible accuracy (through an adaptive refinement). 

c) The solver uses the most simple input. 

d) The solver combines the advances in the mathematics and computer science 
including parallel computations. 

The solver FEARS (Finite Element Adaptive Research Solver), its mathematical 

version FM developed for Univac Series 1100 implements some of the points mentioned 
* 

above . The detailed description of FEARS and the experience with it will be 
published elsewhere. For some information about FEARS and its applications, 
we refer lo | ;> ] , f 3 ] , | 4 ] , [ 5 ]. 


* 


For analysis of the parallelity we reter to [ 1 ]. 
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One of the main aspects of FEARS is the theory of the a-posteriori estimates. 

Some aspects of the a-posteriori error estimates and optimality of the meshes 

were investigated in [6 1 , [ 7 1 , 1 8 ] . 

The error norm ||e|| is approximated by the computable estimator e 

computed through error indicators n^(A) associated to every element A and 

computable locally by knowledge of the finite element solution at the particular 

£ 

element A and its direct neighbors. The effectivitv index 0 = || —| j- expresses 
the quality of the estimator and 0 should be close to one when the error is 
sufficiently small (e.g. 5/1). It is desirable that the estimator z has the 
following two properties: 

(1.1) 0 C L 1 e 1 c u < ” 

with C and C., independent of the solution and the meshes under very general 
conditions. 

(1.2) 0 -> 1 as | |e| | -*• 0 

provided that some additional assumptions about smoothness are made. The present 
paper develops the theory of the estimator which satisfies (1.1) and (1.2), and 
is implemented in FEARS. The energy error norm is assumed and model elasticity 
problem is considered. 

Section 2 consists of some preliminary notions. 

Sections 3 and 4 elaborate on the type of meshes which are adaptively 
const ructed. 

Section 5 deals with the approximation properties of the elements on the 
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admissible meshes. 

Section 6 formulates the model problem (elasticity problem). 

Section 7 develops the estimator and proves (1.1). 

Section 8 proves that the estimator is asymptotically correct, i.e. 

0 -> 1 . 

Section 9 deals with two computational examples and discusses the effectivity 
of the approach. 

The adaptive construction of the meshes is based on the equilibration of 
the error indicators. This principle was theoretically analyzed in [ 7 ] for 
one dimensional problems and its theoretical investigation in the context of 
FKARS will appear elsewhere. 


I 

! 
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2. BASIC NOTATION 


2 

Throughout this article we denote by R the two dimensional real Euclidian 

2 22 1/2 

space with x? , ||x|| = max( | x^| , | Xj | )»| | x | | E = [x^+x^] . Let 

2 

QCR be a bounded set and SQ its boundary. We define 


diam Q = sup ||x-y|| 
x,v6Q 


dist(x,Q) = inf||x-y|| 
y€Q 


and for 


2 

Q i € R 


i 


1,2 


dist(Q ,Q„) = inf | | x-y | | 
x€Q 2 


E 


An index E will denote that the norm ||•| j is used instead of ||‘|| 
, dist (x,Q) = inf||x-v||_ 

y^Q E 

For a p > 0 , Q C is the p-neighborhood of Q 


0 


P 


{x€TR 2 


dist(x,0) 


< P 


} . 


The closure of Q in IR 2 is denoted by Q , int 0 means as usual the interior 
of Q . 

2 

By 7. we denote the set of all two dimensional integers k i (k^,k 2 > , 
kj, i = 1,2 integral. 

2 

Suppose 6 > 0 is a positive real number, then we will write for any k€Z 
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Q k = {x€IR 2 Ik.0 < x. < (k.+l)6 , i = 1,2} 

e 1 i - i - i 

2 

Assume that Z C Z^ is a finite set. Then we denote 
o 



int [ 


U Qgl 

k£Z 9 
o 


We shall assume that Z q is such that 0 is a Lipschitz domain. For brevity, 

o’ 

whenever it cannot lead to misunderstanding we shall write instead of q . 

2 ° 

When we talk of a square S in IR , we shall always suppose that it is closed 

and that its sides are parallel to the coordinate axes, i.e. S is of the form 

[a,a+dl x [b,b+d] for a,b,d£IR, d > 0 . 

As usual, let L ? (Q) = H°(f2) be the space of all square integrable functions 

on SI with the inner product 


(u,v) 


L 2 (fi) 


uvdx, dx = dx^dx 2 




and the corresponding r.orm j 1 
the usual Sobolev space with the norm 


l 2 (w) 


By H (W) , k > 0 integral we denote 


H (Q) 0< u|<k 


l|D a u|l 


L 2 (n) * 


where x = (a^o^), a, > 0 , |a| = 


and 


D° = 


a l a 2 
9x^ 3x z 
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Obviously we have H°(L) = . We will also use the notation j 

! 

S 

I 

| 

i 

f 

in the usual (distributive) way ' 

be the completion of the set of 

We will also deal with functions defined on one dimensional manifolds, 
more precisely on the boundary 3Q or a part T of it. The notation 

L ? (F) = H°(r) has then the obvious meaning. 

m k (i) 

Let F = .U r. , where each T. is a closed side of some 0„ C fl , 

1=1 i 1 < e 

with ILCas? , (i = l,...,m); then we shall write 
H 1,r (il) = {uCH 1 ^.!) u = 0 on IM . 

Obviously (Li) -= H 1,r («) when r = 3S2 and H 1 (fi) = H 1,F (fi) when T = 0 . 

Finally by C°(S2) we denote the space of all continuous functions on Si 
and let 

| | u | | = sup|u(x)I 

C°(n) x€Si 

1 °i 

We will deal later with extensions of functions in H (S'!) , H (SI) and 
H^* ‘ (/) from into a neighborhood of Si . 

Theorem 2.1. T here exists an operator T mapping (Si) into H'*'(S! P ) > 

(where l! 1 is a p-neighborhood of Si ) such that 


l"l\ -l IlnMI , 2 

H (Si) |a[=k J 2 (i) 


We define the support of a function u€L 2 (Si) 

°k k 

and denote it by supp u . Let H (sl)cH (Si) 


all functions having compact support in S2 










i) for a 2 < £ < °°, 0 < ct^ < p and any square SCfi and any 

0 < a < <x we have 
o 


ITull I a 5 CI IuI| 

H 1 (S Jl ) H 1 (S 6 a nn) 


with C independent of a, S, u 

ii) If x€!2 P , x^ft , T ^ ,<ft 


and dist E (x,r) £ distg(x,9ft-f) then 


Tu = 0 on (x)-ft, y = dist_(x,r) 

Z L 


with S (x) being the square with the center in x and lengthside y 


then Tu = 0 on ft P - ft 





iii) If I 

= 

9ft i.e 


Proof. If is enough to prove the theorem in the neighborhood of the boundary 

(i.e., endpoints) of F . In the neighborhood of all other points x€9ft we 

use the classical extension theorem when x f^P and we extend by zero when 

x €T and apply the standard argument with partition of unity. 

The endpoint of T can be located in a vertex of 9ft with the internal 
3 1 

angle y " or t or it can be on straight part of the boundary. 

We will deal only with the case of the vertex in the coordinates origin 

3 

and the internal angle on ^ x . The other cases are analogous. 

Let (see Fig. 2.1) 


ft P = (x€ft f | x^ft, dist (x,9ft-D _< 4 distg(x, T) 1 














<a 
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Obviously ^*** consists of a sector with the lines and as 

its boundary. Let V be a symmetric (with respect to the origin) sector 

% 'V, % 

with the boundary S 1 and and let be the line symmetric to 

. % 

(with respect to <10 - r) . Finally let VCS! (respectively VCO) be the sector 

bounded by 10 - T and S* (respectively f>^) and W be the sector bounded 

^1 

by S and F . 

1 r 

Assume now that u €H ’ ( 12 ) . By an afine transformation, we construct 
w£H (V) such that w = u on - T , w(x) = u(x) for x€S^, x€S^, and 

|jx||p = ||x||g . Obviously 


($) 


i x (h 


The extention Tu on be now the reflection of w (with respect to 

>' -f) . It is readilv seen that 


i I Tu , 


V(* 


ibv) 


1 1 ^ 
H J (V) 


and for 

x €S 1 , x 6.^, I U! | £ = I I x I i £ 

we have 


(Tu) (x) = ti(x) 


We extend u on by zero. 

Bv an afine transformation of W on V we can easily continue a function 
v such that v = u on and v = 0 on !j_ and 
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H 1 (V) 


! u i 


H l (W) 


Now let Tu on 1 be the symmetric image of v . It is easy to see 

that our construction has all properties of the extension formulated in the 
theorem when is chosen sutticientlv large. 


► 

: j 

i I 
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3. THE ME SH AND ITS BASIC PROPERTIES 

Wo will introduce now a class of partitions of SI 

'o ’ 

We define a mesh V(i',) = {A } as a finite collection of closed squares 
A*Cii of various sizes with sides parallel to the coordinate axes, and which 
are generated by the following recursive rules. 

k 

i) The squares (.Q^l > k€ Z , create a mesh. 

ii) If {A 1 } , i =],...,m is a mesh, then a new mesh is obtained if any 
A 1 is subdivided into four congruent squares of half the side length of A* . 

Any & X €D will be called an element and its sides the edges. The 
vertices of the elements will be called the nodes. A node P will be called a 
regular node if either P6 3.Q or P is a vertex of four different elements. 
Otherwise P is an i rregular mode. By P(P) we denote the set of all nodes 
of V and by R (P)CP(P) the set of all regular nodes. Finally let h(P) = 

ma^ diam A . 

A€P 

Figure 3.1 shows an example of a mesh. The irregular nodes are marked by 
a cross. 



Figure 3.1. An example of a mesh. 
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Let D ’.A 1 ) , i = Then obviously Jj A 1 = Si . We shall 

i = l 

denote by AI(D) , the subspace of all continuous functions on fi which are 
Individually bilinear on each A* , i = 1,...,m . It is clear that 

M(P)CH l LA . 

We will always assume that at least four different £\€V lie in every 

Q- • 

Now we will analyze further the basic properties of the meshes introduced 
above. 

LEMMA 3.1 . Let A ' , A" € P 

1) Assume t hat £ 1 f) A" 1 0 and diam A' <- diam A" . Then one and only 

one o f tiie following stat em ents holds. 

i) 

i i) A 1 D " is just one poin t being a common vertex of A* and A" . 

iii) A' fl.A" is__an edge of A' and is contained in an edge. T say. 

of . If x', x" are any two endpoints of (A'PlA") and F respectively, 

then !x'-x"| is an integral multiple of diam A' . 

2) If P €P(D) , P ^R(V ) (i.o., P is an irregular mode) and P is a 

vertex of A' € V , then there exists A"C V such that 

il P € A" 

Li) diam A" diam A' 
iii) A” ft, A’ is an edge of . 

3) diamA"/diamA' = 2 s with s an integer. 

The lemma can he easily proven by induction. 

LEMMA 3.2. Suppose t hat AjG’jC*_._Th en at least one vertex of A is a 









regular node and If an edge of A is contained in 30 for some k€Z , 

___________ 'Q _ o 

k 

then at least one vertex of A which is not on is a regular node. 

Proof. Observe that any regular node always remains regular when our 
recurrent construction is implemented. The lemma follows now easily since at 
each step the midpoint of the subdivided element becomes a regular node and it is 
a vertex of all four new elements created at that step. 

As seen in Figure 3.1 there could be an element A€V , such that only one 
of its four vertices is a regular node. 

LEMMA 3.3. Let V be a mesh and P€R(P) . Then there exists 
Vp€M(P) such that v^ (P) = 1 and v^CQ) = 0 for an y Q€R(P) , Q 4 P . 


Proof. We first note that it suffices to define v at the vertex of each 
- P 

&€V . Let V = {A 1 } , i = l,...,m . Assume that we have denumerated the 
elements so that diam A 1 >_ diam A 1+ ^ i = l,...,m-l. We prove now our lemma 
by induction with respect to i . 

First let us observe that all four vertices of A 1 are regular nodes. 

If one were irregular then by lemma 2.1, there exists A€V such that 

% 1 1 
diam A > diam A . This is a contradiction because A is the largest element. 

Suppose now that j = 1 , or j > 1 and v has already been constructed 
i=j-l . . P 

1 1 T 

on U A Consider now the vertices of A . If Q is a vertex of A J 

i=l 

and 

i) Q is a regular node, we define v^CQ) = ^ for Q = P and v^CQ) = 0 
for Q / P . 

ii) Q is an irregular node, then by Lemma 2.1, Q€A' where 

diam A' "* diam A-' . By induction assumption v^ was already defined 

on A' and so v (Q) is defined too. We therefore construct v 

i = l K 

on A^ with desired property on U A Let us remark that 

i=l 
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i=j-l . 

if 0 £ U A then the value v (0) = 1 respect ively 0 is the same as in 
i = i p ‘ 

the previous phase. 

LEMMA 3.4. Let D he a mesh and u£M(P) be sucti that u(Q) = 0 for any 
Q€ R(P) then_ ju = ()_ . 

Proof. Let {A*j be numerated as in the proof of Lemma 3.3. The lemma 

will be proven by induction with respect to i . Because all vertices of A^ 

are regular nodes we have u = 0 on A^. Let now i = l,...,j-l < m and consider 

u on A - * . If the vertex P of is a regular node then by assumption 

u(Pf = I) . If P is irregular node then by iemma 2.1, we have also P f i' 

with diam t\ dinm A" . i.e.. A* = , for some k '■ j . By induction 

i = .i-l j 

u = 0 on .U, A and therefore u(P) = 0 . So u = 0 in all vertices 
of and so u = 0 on A * and lemma is proven. 

Lemmas 3.3 and 3.4 show that the function u€M(P) is uniquely defined 
hv its values at the regular nodal points. 

Lemma I. t and arguments analogous to those used in the proof of Lemma 
3.3 yield 

LEMMA 1.3. Let l’ he a mesh. Then 

i) The set of functions 1 v |P €R(V)} creates a basis for M(V) 

i i ) v ■ D 
P 

i i i) ) v = 1 

P € R(V) r 

Definition 3.J: The set = supn v , J’€R(P ) will be called a star 

: v:—~r—77: _JP_P___ 

associated to the node P or briefly a P-star. 

_ f 1 

P 


Lemma 3.3 (iii) vie his readily that U 

P €R(V) 
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LEMMA 3.6 . Let t,€V , P€R(£>) and A Hint m p 4 0 . Then iCin p . 


Proof. Assume on the contrary that A<^iu . Then there exists an open 

set SCA , Sfl 'ii = 0 . Because v £M(D) , v is bilinear on A . But 
P P P 

v =0 on S , and hence v =0 on A . This is a contradiction because 
P P 

we assumed that A fl int w ^ 0 . 

P 

LEMMA 3.7. The set o is connected in the sense that for anv two elements 


A’, A"C ui there exists a sequence of elements A 1 = A„,A.,...A = A" such 

n M 0 1 n 


that 


i) A.Cio 
-—1—_P- 

ii) ^ A. +1 ^ 0 and is t he edge either A^ or A^ + ^ 


S i j _ . 

Pr oof . Let = (J A J , A 1 € y . Assume that we have enumerated the 


j-l 


'Kwnf.s A J of 


so that diam A'' > diam A^ + 1 , ) - l,...,s-j . It is 


’ 1 

obvious that we can restrict ourself to the case when A = A 

First we prove that one of the vertices of A^ is the node P . Suppose 
that the node P is not a vertex of A^ . Let Q be any vertex of A^ . If 


Q is regular, then v (Q) = 0 . If Q is irregular, then there exists 

f\*€V with Q € A* and diam A* > diam A^ . Thus A* flint w = i)> and so v = 0 

P P 

on A* by lemma 3.6 and it follows that V p (Q) = 0 , so v p = 0 on A^ 

which is a contradiction. Now we prove the lemma by introduction with respect 

k 1 

to j . Assume therefore that we are able to connect A with A and 

k+1 k+1 

consider the element A . If a vertex of A is the node P , then 

k+1 1 

A can be connected with A because n is a Lipschitz domain. If all 

k+1 k+lj. 

vertices of A were regular nodes different from P then A <tw by 

’ P 

lemma 3.6. 


So we need onlv consider the case where a vertex R of A 


k+1 


is irregular 






and v (R) ^ 0 . Bv lemma 3.1 there exists A* € V such that diam A* > 

P 

k+1 _ k+1 W+1 

diam A , A* HA is an edge of A and R € A* . Thus A*Cu and 

P 

so A* = a' , j ' k . So A^ + ^ can he connected with A^ and therefore 

with A^ and lemma is proven. 

Lemma 3.7 shows that int w is a domain. 

P 

So far we have not made any restrictions concerning the mesh V . In 
the next section we will analyze the family of K-meshes, which play an 
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A. THE K-MESH 


Del in it ion A.I. Lot K • 0 , real, A mesh V will be called a K-mesh 


if for any I 1 € R(V) 


(A.1) diam , < K inf diam A' 

P A'€V 

A' C i. 


The definition has a clear sense because of Lemma 3.6. 

We conjecture that definition A.l is equivalent to 3K* > 0 such that 

j~ -i i n diam A , 

for all A € V , sup —--—< K* 

. n . i . , c , diam A — 

Alla is an edge of A 

Everywhere in what will follow we shall assume that we deal only with 
K-meshes. We mostly will not mention it explicitly. 

LE MMA A.l. Suppose V is a K-mesh. Then there exists numbers M r N 
depending only on K such that 


O_I f_ P€R(P) then the P-star consists at most of N differen t 

elements of V . 

ii) If A '€V then A'Cm^ for at most M different F€R(V) . 

iii) I f_ then there are at most AK + A elements A"CP such 

that A ' PI A" i 0 . 


Proof. 

i) The star u> can be contained in a square S with its side diam u > 

P P P 

So for the number N of elements contained in oi we have the simple estimate 

P P 


area m fdiam m 1 

__p __ p 

inf (area A') - [inf (diam A ')]^~ — 

A '£ -o A 
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The node F must he a vertex of some A"C w satisfying (4.2). The number 

of such elements is obviously bounded by 


area Q 

_L_ 

.di am A' 

' V ' 


4K 2 (K-l/2) 2 


Because not more than 4 regular nodes could be on any element, we see that 


M < 16R 2 (K-1/2)" 


iii) Any vertex of A’ can be a vertex of at most three other elements 

and therefore by lemma 3.1 it is sufficient to bound the number of elements 

A" with diam A" _< diam A’ such that A^flA’ is an edge of A" contained 

in some edge of A' , say P . Fix this T , and suppose that there are q 

such \" . q is finite since by the definition of V , there is only a 

finite' number of elements. Then for at least one such A" , A say , diam 

— diam A' . Since y v = 1 by lemma 3.5, there exists P€R(V) such 
q P€R(D) p 


that v is not identically zero on a Da 1 

p 

A ' C u) . But 


and therefore both A and 


q diam A < diam A' < diam m < K diam A 

p _ 


Hence 


q ; K 

and the lemma easily fellows. 

LEMMA 4.2. There exists a number L (depending only on K) such that 


o><? 
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R(V) can he partitioned into j < 1. sets x j > j = such that if 


I\o€ , 


J 


5 s Q then int 0 int Wq = 0 


Proof. Suppose P 4 Q , P,o€R(V) and int w^PI int Wq 4 0 • By lemma 
3.6 int l> flint au must contain the entire interior of at least one element. 

P Q 

Bv Lemma 4.1 ,o contains at most N elements and each of these elements 
P 

can be contained in o , for at most M-l nodes P' € R(V) , P 4 P' • There- 

P 

fore there can be at most N(M-l) regular nodes Q such that int LOpflint 

We shall construct now the sets Xj by the following recursive procedure. 

Let P , P ? ,..,P r be some enumeration of the regular nodes. Set = {P^K 

l-l l-l 

Suppose that we have already defined sets v^ ,...,Xg for some >_ 1 , 


l-l 


l ■ H. 1. [f for P^€R(P) and some 1 £ k £ S^_^ , int u> flint - 0 


for 


all Q€ 


l-l 


, then choose k to be minimal, and set 


/ = J for t €(],... ,S^ - (k) 


4 - 4' lui, v 


and define = S^_^ . Otherwise set 


l £-1 

, t = for 1 < t < S 


l-l 


>S £-1 + 1 = {P £ } 


w i tii = S^_ j + 1 . 

Now from the first part of our proof we see that N(M-l) + 1 and so 

L < N(>1— 1 ) + 1 and the lemma is proven. 
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Remark. As we mentioned in the introduction, this paper develops the 
basic ideas of [1] . The lemma 4.1 relates to the intersection index and the 
lemma 4.2 to the overlap index as introduced there. 

Given an element A€P we will always enumerate its vertices as shown 
in Figure 4.2. 



Figure 4.2. The Numeration of the Vertices 


Let P€R(D) and w its star. The node P is a vertex of at most 4 

P 

elements. We will denote by A cm the element for which the vertex number 

r r 

given by Fig. 4.2 has the minimal value. This rule associates to every P a 
unique element. 

Let now m be a P-star. Then by J we denote the invertible affine 
P P 

transformation taking P to the origin and (A^> = [0,1] * [1,0] . Further 

let $ = J (to ) . 4> will be called a standard P-star and we shall call 

P P P P 

members of S = {J (A)jA€t? , AC“ } the standard P elements of 4> and 
P P P P 

if no confusion arises denote them also by A . Note that <S> can equal 4> 

P Q 

for P 4 Q and yet 5 y S . 

P Q 

LEMMA 4.3 . 

i) int i> is a domain. 

_E_ 

ii) There are not more than Z = A(K) possibilities for S as P ranges 
- p- 



over R(V) . 
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Proo f . 


i) Follows Immediately hv lemma 3.7, 


ii) The result will follow if we can show that there is a finite collection 

of squares S,,...,S ... such that {J (A) |a€D AC“ } = {S.....S \ } 

1 -n(k) p 1 p 1 ir(k) 

To this end we show by induction that the vertices of each J (A) must 

. _ P 
-([^.2 kl int +l) 

have coordinates of the form (kp,^n) where n = 2 v s 2 J k,£ integral 

|k|,j£| £ 2K^ and that diam J^CA) = 2^n for some T = 1,2,3,... Here [']^ nt 

denotes the integral part. 

Now if A , ACw we can construct a sequence 


A = A ,. . , A = A 
p o’ n 


having the properties mentioned in Lemma 3.7. 

Clearly any vertex of J(A^) satisfies the above inductive hypothesis. 

So suppose it holds for J(A^),...,J(A^) 0 £ i ^ n-1 . Then by Lemma 3.1 

(3) diam = 2 s diam A^ , s integral, giving diam = 2 S diam J^CA^) = 

2^ +S n for some integral A,s . Since A..,C cu we have as in 4.2 diam A.. n > 

l+l p i+1 — 

diam A /k . This gives diam(J (A...)) > — and so K + s > 0 and so we 
p p l+l — K 

conclude that diam J^(A i+ ^) = 2% with t >_ 1 integral. 

Suppose for the moment that A^ + ^ 2L diam A^ . Appealing to Lemma 3.1 (1) 
we see that the two vertices of J (A, which are the end points of 

p l+i 

J (4ui>nj (A.) must satisfy the induction hypotheses. Since J (A,...) 

P I' 1 p 1 p It1 

is a square, it follows also for the other two vertices. The bound on 
|k|,|£| is a consequence of diam <t £ 2K . The case diam A^ <_ diam A^ + j 
follows by a similar argument. 


COROLLARY: There is not more than Z(K) of different possible domains Q 


f 
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Given t>€V we denote 

(4.3) Q*(A) = U{u lp |p€R(P) , ACa) p } 

As in the proof of Lemma 4.1 we have Q*(i)£Q(A), where Q(A) is a square 
with the center at the middle of A and 

(4.4) diam Q(A) 2(K-l/2) diam A 

LEMMA 4.4 

i) If A',A"€P and A',A"CQ*(A) then 


_1 

(4.5) _ diam A ' > diam A 

and 

(4.6) diam A' > —If diam A" 

___ _ K~ _ 

ii) Int Q*(A) is a domain . 

iii) If P€R(V) and w p CQ*(A) then on A 

(4.7) |D^v j < K(diam A) ^ . 

P - 
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Proof. 


hence 


Let A'C Q*(A) 


then for some P' €R(P) 


we have 


A,A'€ «•> , 


(4.9) 


diam A' 


> — dlam w , 
- K p' 


>_ — diam A 


Further we have for A"coj „ 

P 


(4.10) diam A > v diam w ,, > -i diam A" 
— K. p — K. 


(4.9) and combination of (4.9) and (4.10) yield (4.5) and (4.6). 

ii) We have to prove only that Q*(A) is connected. This follows 
immediately from the Lemma 3.7. 


iii) From Lemma 3.5 (ii) and (iii) we have 0 _< v <_ 1 and so obviously 


and 


!l) 1 v 


[ min diam A] ^ 

AC u> 

P 



[ min diam A 

AC to 
P 


n 


(4.9)yields the Lemma. 

s o 

Let K < 2 . Denote 


Q**(A) 


(x €Q*(A)|dist(x,3Q*-3Q) 


diam A 

> 

2 ° 


} 







LEMMA 4.5. 


i) Let 'o CQ*(A) . Then P€Q**(A) . 
___ P_ __ 

ii) ACQ** (A) 

iii) int Q**(A) is a domain. 


Proof. 


s„+2 


i) Assume that P^Q**(A) then dist (P,3Q*-3£i) < diam A/2 ° 

Because v^ = 0 on 3Q*-3Q , (4.7) leads obviously to a contradiction. 

ii) Assume that A<^Q**(A) . Then there is a vertex P* of A such 

that dist (P*,3Q*-9f2) < (diam A) 2 ^ So+ ^) . Obviously w = £ v = 1 

WjPQ*(a) P 

on A . By the same argument as leading to (4.7) we see that |D^w| <_ K(diam A) 

Because w = 0 on 3Q*-3fi we have a contradiction. 

iii) Let xf Q**(A) , x€A C oj . Then by Lemma 3.7 there exists sequence 

o p 

n 

A A = A such that V =int U A. is a domain. Because diam A. > 

° 1 i-0 1 J ■ 

— diam A by (4.2) it is readily seen that int Q**DV is a domain. This leads 
K 

immediately to the desired result. 

Let if* = Jp(Q*(A)) (analogously as <J> p ) and if** = J p (Q**(A)). Then we have 


LEMM A 4.6. There exis ts not more than Z*(k) (respectively Z**(k)) 
possible domains int if* (respectively int if**) . 


Proof. The first part of the lemma follows from Lemma 4.3 and is a 
corallary. The second part follows by analogous argument when realizing that 
if** is composed by squares with diam n/4 when n was introduced in the proof 


of Lemma 4.3. 
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5. THE APPROXIMATION PROPERTIES OF M(P) 


Ln this section we will analyze the approximation properties of M(D) . 
Let R be the mapping of C°(0) onto M(t?) such that (Itu) (P) - u(P) , 
p€R(P) . By Lemma 3.5 we can write 

n u = 7 u(p)v 

P€R(V) P 

Further define the operator ^ ^ mapping C (fi) into M(P) by 

(5.1) n Q * (A) (u) = )' u(P)v 

P€R(0) P 

io CQ*( '.) 

Clearly supp(li^ ^^(u))C0*(A) and 

(5.2) ri Q * (A) (u) = P(u) on A . 


For given A€0 we define as the invertible affine transformation 

2 2 

of TR onto IR taking A into standard unit square [0,1] * [0,1] = S . 

LEMMA 5.1. There exists a constant C dependent only on K such that 
for any w€H^(n) and £=0,1 


(5.3) 


I woJ. 


(n Q * (A W 


I £ 

H 4 (S) 


< CI IwoJ 


-li 


H^(J a Q**(A)) 
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Proof. By Sobolev imbedding theorem we have 


(5.4) 


| woJ 


- 1 , 


c°(J a q**(a>) 


< CIIwoJ 


- 1 , 


H 2 (J a Q**(A)) 


Because of Lemma 4.6 there are at most different domains J A Q**(A) , 
where Z**(K) depends on K only. This shows that C in (5.4) depends on 
K only. 

Using (5.1), Lemma 4.4 and Lemma 4.5 i) we get 


(5.5) 


(n 


Q*(A) 


w)°j" 


I l 
n (S ) 


< CI IwoJ 


-li 


H 2 (J A Q**(A)) 


where C depends on K only. 

Since by Lemma 4.5 ACQ**(A) 


(5.6) 



H^(s) 


< 


llwoj" 1 !! 2 

h^(J a q**(a)) 


Combining (5.5) and (5.6) we get the Lemma. 


LEMMA 5.2. There exists a constant C (dependent only on K) such that 


for any w€H (fi) and £ = 0,1 


(5.7) IwoJ" 1 - (II Q * (A) w)oJ iClwor 1 ) 

H^(s) A H Z (J A Q**(A)) 

Proof. Suppose z is any function bilinear on Q** . Then on A we 
Q*(A) 

have IP z = z and therefore 
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(w+zOuJ. 1 - (!1 Q * (A) (w+z))°J 1 = woj 1 
A A A 


(ri^W 1 

A 


and hence using lemira 5.1 we get 


Iwoj" 1 - (n Q * (A w)°J 1 


■1 .i 


a w/ -'A 1 1 Q - ^ * n H l Wo -^A “ z l I 2 

A H V '(s) z A H(J Q**(A)) 

A 


where z is an arbitrary bilinear function on J^Q**(A) 


Because by lemma 4.6 there is only a finite number (depending on K only) 
of different J^Q**(A) , and by lemma 4.5 any int J^Q**(A) is a domain, we 
have (see e.g. ( 2 1, pp. ) that 


Inf[(w°J *-z|| c|w<tJ X | 

H / (J A Q**(A)) H^(J a Q**(A)) 

The theorem follows immediately. 

Now we have 

THEOREM 5.3. There exists a constant C depending only on K such that 
for any A CP, uCH 2 (i2) and l = 0,1 

(5.8) I 1 11 -rr u | I < C(diam A) 2 ? | | u | | 

H '(A) H (Q**(A)) 


Proof . Using (5.2) we obtain (5.8) from lemma 5.2 by standard scaling 
argument. 


THEOREM 5.4. Let uCH 2 («) then for H = 0,1 






with C depending only on K . 


Proof. We have 

||u-nu || 2 = I i!u-n u |! 2 £c£|MI 2 i h 2(2_fc) (P)< 

H (ft) A €0 H*(A) A€D H tQ**(A)) 

1 C J ||u|| 2 h 2{2 ~ n (V) < C l X(A) | |u| | 2 h 2(2 "°(P) 

h€V H Z (Q*(A)) h€V H*(A) 

where A(A) is the number of Q*(A’) such that ACQ*(A') . Using lemma 4.1 
we see that A(A) £ MN £ C(k) and so theorem 5.4 is proven. 

2 

Remark: In Theorems 5.3 and 5.4, the restriction u^H (ft) can be weakened 
to u€H 2 (A) for every A€P . 

Theorem 5.4 shows that M (V) has the same basic "interpolation" properties 
as the usual finite element spaces. The spaces M(P) are more flexible than the 
usual spaces defined on quadrilateral meshes. The space M(D) allows us to 
make a refinement and still keep square elements. The restriction to K-meshes 
is from a practical point not important. A more essential restriction is that 
we deal only with squares. How to overcome this resctiction with respect to the 
implementation and the theory will be discussed elsewhere. 

The use of the spaces M(V) does seem to have a major advantage over tri¬ 
angular elements because of programming and data management simplicity, expecially 
when some form of automatic or adaptive mesh generator is envisaged. One mani¬ 
festation of this is that the element can be uniquely defined by the binary 
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expansion of the coordinates of its center, the length of such an expansion 
indicating the size of the element. Since all the elements have moduli, a 
scaling factor the same geometric shape, the calculation of the stiffness 
matrix, etc. is simplified. In addition this seems to be important for the 
practical effectiveness of the error estimation. 

We shall analyze now the approximation properties of M(P) when u€H^(Q). 

2 

Before being that, we introduce some notations. By p(x), x €1R we denote 
a molifier, a function with all derivatives continuous p(x) >_ 0 , p(x) = 0 
for |xl > 1 , u(0) = 1 , and 


p(x)dx = A . For e > 0 , let y^(x) 


IR, 


c 2 A 


d(x/e) 


Let ftCIR 2 be an arbitrary bounded domain and ff its p-neighborhood. 


For u defined on 0 , we put 


u = u*p , c < p 
c e 


Obviously u is defined on fi . 

J e 

2 

Further for any t €fR , |t| <_ 1, e < p let 


u te (x) = u(x+tc) 


Then u tC is also defined on , and we have 


(5.10) ||u || 2 <Cc 1 ||u|j 

E h (n) h 1 (fi p ) 


l u f - u l i (, 1 Cf;1 ? 'l l u l I , _ » O' = 0,1 

H (0) H (if) 
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Iu t£ -a|| 1 Ce 1_2 ||u|| 

H (fi) H (ST) 


with C independent on f! , u and e 


Let us prove now 

LEMMA 5.5. For every ?€R(V) , let a function w 6H 2 (fi) be given. Let 
- — - - - ___E- 

A € V and 4Cw . If 
Po 


w = n y v w 

P€R(0) p p 


then for l = 0,1 


(5.11) 



/(A) 


C[(diam A) 2 ^ 


w p M 2 

P o H (Q**(A)> 


+ / y | |--w +w | | . (diam A) * 

P€R(P)i=o p o p hVU)) 

Ae. 

P 


Propf . On A we have 

w » n y vw + n y v[-w+v] 

P€R(5) p p o p€R(a) p P o p 

= n Q * (A) w + n Q * (A ' l v (-w +w ] 
p o P€RiA) p P o P 


Using lemma 4.4 and Leibnitz's rule, we have for l = 0,1,2 
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(3.12) 


l I 

A ' n Q**(A)^0 




(q**(a) n a' ) 


C ) ||-w +w |j . (diam A) 

i=0 p o P H 1 (Q**(A)) 


i-P 


Therefore on A 


„ - ,,Q*(A) 

—w = w -n w 


7 V [-W +w ] 

P€R(A) P p o P 


[n Q ^ A) ( 


7 v f-w +w ] 
P€R(A) p p o P 


7 V [ -W +W I ) 

P€R(A) P P o P 


Applying theorem 5.3 and noting the remark after theorem 5.4, we get the lemma. 
LEMMA 5.6. Let v€H 1 (IR 2 ). Associate to every P€R(V) a sector 

9 

t € IR , |t I r 1 and a number A such that 

P P_—r __ 


>, 

P 


min(diam A )' 

A€ - 

P 


1 

K.512 


Let further 


t 3 
, P P, 


1 

-J?. 

64 


and 


w = n 7 v w 

P€R(V) p p 
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Then for any A€P and i, = 0,1 


|v-W| I <_ C(diam A) 1 £ | |v| | . „ 

H (A) H (Q**(A)) P 


(5.13) p = max A 


u» p CQ*(A) 


Proof. Let ACw and ACw . Then using lemma 4.4 we have 

- P P 

r o 

A p ,X p £ C(K) diam A. Using (5.10) we get for i = 0,1,2 


w -w . 

P p 1 ' i 


< C(diam A) 1 i I|v| 


o K (Q**(A)) 


H 1 ([Q**(A)] 2p ) 


By lemma 5.5 we get readily the lemma. 


LEMMA 5.7 . Let p be defined by (5.13), then 


[q**(a)] 4d n«cq*(A) 


The lemma follows easily from the definition. 

THEOREM 5.8 . Let u€H^(fl) respectively H^ ,r (fl) . Then for any V 

1 r 

with h(0) _< there exists w€M(V) respectively M(P)Dh ’ (f2) such that 
f or all h€V and l = 0,1 

(5.14) [ |ii—w| | < C(diam A ) 1 | | u | | 

H^(A) H (Q*(A)) 



J 
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with C independent of V, u and A . 


Proof. The extension of u onto a neighborhood ft 0 , has the properties 

2p o 

listed in theorem 2.1. Denote by [Q**(A)] the 2P o -neighborhood of Q**(A) 

p = min[o,p/2Sl where p is defined bv (5.13) and 8 was introduced in 
o 

Theorem 2.1. Then bv lemma 5.7 and theorem 2.1 


IN! x 2 p - I i 4p I C M U II i 

U (fQ**(A) ] p ) H i <[Q**(A)] P Dft) 1>[Q*(A) ] 

Let u€H^(ft). Select t = 0 in Lemma 5.6 for all P€R(V) . Then 

P 

for W constructed in Lemma 5.6 we have W€M(0) and 


|u~W|| < C||u| 

h l (a) 


! 1 2p i c 11 u 11 L 

H i [(Q**(A) ZP ] H i [Q*(A)] 


and the first part of the theorem is proven. 

1 F 

We have now to show that when u€H ’ (ft) then we can choose W so 

that W(P) = 0 for all P€R(D), P€T . To every P€R(D), P^F , we take 

t = 0 . Tf P €. 3ft, P € F , p€?V then we take t the outward unit normal 
P P 

(if P is a corner of ft then the normal is bisecting the outside angle). 

If P€lf’ then it is easy to see that we can select vector t (pointing in 

to the sector ft P y . in Fig. 2.1) and not necessarily of the unit length such 

that w (P) = 0 also. This finishes the proof. 

P 

Finally we prove 

i ° I i r 

THEOREM 5.9. Let u€H (ft) [respectively H (ft), H ’ (ft))]. Then there 


exists w€M(P) respectively M(P)0H^tR) respectively M(P)nH 1,r (ft) such 


that 
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(5.15) l |!v (o-w)|| 2 £ C||u|| 2 

P€R{V) p _ H (fl) _ H (Q) 

with the constant C dependent only on K . 

Proof . By Theorem 5.8 we see that for any A€P and i = 0,1 

(5.16) ||u-w|| £ C(diam) 1 ^||u|j 

H (A) H (Q*(A)) 

and applying Leibnitz's rule and lemma 4.4, we obtain 

(5.17) ||v (u- W )|| iC||u|| 

P H i (A) n (Q*(A)) 

for all P€R(V ) such that w^CQ*(A) . By lemma 4.1 there are not more than 
M(K) such nodes so we have 


P€R(V) 


l I |v (u-w)I I x 
(V) P H A (A) 


l 

P€R(P) 

“> CQ*(A) 


V (u-w) 
p 


I 1 

h\a) 


< C(K)I IuI 


H (Q*(A)) 


and therefore 

(5.18) l | |v (u-w)|| 2 < C(K) l ||u|| 2 

p€R(P) p H («) H (Q*(A)) 

By lemma 4.1 there are not more than C(K) different A'£0 such that 
Q*(A')3 A and so (5.18) yields (5.15). 
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Remark to Theorem 5.9. Assume now that int supp u = Q* . Then 


supp w = U w 

P€R(V) P 

o- np.>M0 

p 


This observation follows immediately from (5.16). 


Remark: Theorem 5.9 is closely related to (3.1) of [ 9 ] which is an 

essential part of the a-posteriori error analysis. 
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6. THE MODEL PROBLEM, ITS FINITE ELEMENT SOLUTION 


AND THE BASIC A-POSTERIORI ESTIMATE 


As a model problem we will discuss the case of plane elasticity for a 

body, homogeneous and isotropic on every Q making up the domain ft . 

0 

Let H 1 (ft)CH.CH 1 (ft) i = 1,2 where H. = H 1,ri (ft) as defined in 
Section 2. 

Let Hq = H^ x H 2 , with u = (u^,u 2 ) and consider on Hq x Hq the 
following bilinear form 


( 6 . 1 ) 


B(u,v) = B(u 1 ,u 9 ;v 1 ,v 2 ) = 


3u 3v 3u„ 3v 

(A+2p) ( -r— 11 + ~ "r —■) + 

3x^ 3Xj Sx^ 9 x 2 


3u^ 3u„ 3v, 3v 

+ U ^3x 2 + 3x^''3 x 2 ' 3x^ 


~ dV dV du 

+ + M-r-i 


3v„ 3u„ 

3 v - 

2 2 

1 

3x 2 3x 2 

3xi 

1 constant 

on ( 


H 


dx 


the constants A , u are the usual Lame constants. 

We will assume that there exist constants 0 < c r c 2 < 00 such that 


for any u € H 


0 


(6.2) C ||u|| 1 B(u,u) 1 C !|u|| 

H 1 (ft) H 1 (ft) 


wi th 


2 

H X (n) 



H X (ft) 


IK 


2 

H 2 (n) 


V 
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ObviousIv B(u,v)‘= B(v,u) and |B(u,v)| < c||u|| ||v|| • There- 

H (SO H (ft) 

fore on H^, B(u,v) is a scalar product with the energy norm 

(6.3) j ||uj j j 2 = B(u,u) 

(6.2) shows that the enorgv norm | | !u| I I" is equivalent with | j u j j . 

H <fi) 

The problem P(H w,g), wCl^O) x H l (S2) , g€ll°(fl) x H°(fl) consists 
of finding u€h'(:..>) x H^Cl) such that u - w € and 

(6.4) B ( u , v) = (g.v) , WvCH^ 

where we have written 


(g,v) 


(g ] v 1 + g 2 V' 2 )dx 


It follows by the standard theory that u exists and is uniquely determined. 

The function u wiLl he called the exact so luti on of the problem P . 

Let now ,M (P) = M(P) OH i - 1,2 and M o (0) = M^P) x MjfD) . 

') 

Assuming that wC(M(P)J^ in the problem P(H ,w,g) , the finite element 
solution U of P with respect to M consists of finding U€M(P) x M(P) 
such that (P) and 


(6.5) B(U,v) * (g,v), \/v€M Q (V) 


Just as for the exact solution it follows that the finite element solution U 
exists and is uniquely determined. 
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Finally we denote by e = U - u the error of the finite element solution 


We will be interested in an a-posteriori estimate of the norm 


|' jej || or some norm equivalent to it. We will design an estimator E 
depending only on the known finite element solution — which will be related 
to the error norm. As estimator E is called an upper respective to lower 
estimato r if there exist constants A^ respectively A^ independent of V 
and u , U such that 


1 M I i \,E 

respective 

A l E < |le|| 


THE O REM 6.1 . Let wCM(P) u be the exact solution of the problem 


P(H o> w,g) and U its finite element solution with respect to M . Then 
there exist strictly positive constants C^, depending only on K , 

u , A and 0 such that 


( 6 . 6 ) 


“0 


JL B(n n* n n) i 


PCR(O) 


< C 1 I B(n ) 

“ 1 P€R(V) p p 


where 


(6.7) 


n €H(i' ) = {v€H n |v=0 on 0-w } 
P P O' p 


and 
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Proof. It follows from the definition of U and u that 


B (e ; v) = 0, Vv^M 


and therefore 


B(e,e) = B(e ,e-v) = B(e, }’ v (e-v)) = l B(n ,v (e-v)) 

P€R(P) p FCR t V ) p p 


Using Schwarz's inequality we get 


B (e , e) 


P€7 1(0) 


I ! !v (e-v)! | I 
P 


t l 

P€R(V) 


2 


1/2 

] 


P€R(V) 


v (e-v 
P 



By (6.2) 


(e-v) ; 1 | “ I. C l I v (e-v) ] \ ~^ 


H (SI) 


and therefore by Theorem S.9 for a proper choice of v 


/. 

pez(D) 


i v (e- v) 

P 


i 2 


Cl iel 


h’oi) 


2 


Hence we have 


P€R(P) 


which proves right hand side of(6.6). 
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Let us prove now the left hand side of (6.6). Define 

w. = y n , i = 1,...,J , where 

J pe \i P 

v. are the sets introduced in lemma 4.2. We have now 
J 

(6.9) B(e,w.) = l B(e,n ) = J B(n ,n) 

J P £, P pg x P P 

j 'j 

Because for P ^ Q , P,Q£X| lint (.supp n ) Dint (supp rig) = 0 we have 

B(n , n,J = 0 3nd hence 
P Q 

(6.10) B(w., w.) = l B(n ,n ) 

1 J PC X J 1 

Further 

|B(C.,W.)| 1 ilk’ll! I I |Wj !! I 

and hence using (6.9) and (6.10) we get 

11 !”,i I ! z - 5 III" 111 2 1 IIMII III-,III 

P€X. P J 

and hence 

in-in 2 :. i!i».in 2 - y 11in jii 2 . 

P€X : 


Because j ranges over 1 I we get 


2 


( 6 . 11 ) 


>: Mini 

P€R (£>) p 


By 


lemma 4.2 -I < L(k) and therefore 
We have proven in Theorem 6 that 


(6.11) proves the left side of 
the estimator 


( 6 . 6 ). 


£ (10 = l B(n , 

P€R(V ) p 


n P ) 


is simultaneously on upper and lower estimators. The individual terms Hp 
are determined locally on separate stars w . Let us underline that al¬ 
though the unknown error e is present in the definition of np in (6.8) 
we have not to know it. In fact 


B(e,v) = B(U-u,v) = B(U,v) - B(u,v) = B(U,v) - (q,v) 

Remark. The proof of the theorem 6.1 follows very closely the ideas 
in [ 9 ] . 

We assumed in Theorem 6.1 that w€M(P) . It is obvious that 
w = (wpiv^) influences the solution only by its values at the boundary 
, more exactly on f . In general when w^M(P) we replace w by 

f\j __ 

w£M(P) and estimate the norm of the solution of the problem P(H o ,w-w,0) . 
Usually it is easy to explicitly construct a function z€ H^(ft) x H^(ft) 
z i = (w-w) j on Fj , i = 1,2 and the desired estimate is then simply ||UIII 
In practical cases we can expect that IIUII I is much smaller than |||e||| . 
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7. THE A-P0STERI0R1 ERROR ESTIMATE 


Let us denote by Q any of the squares comprising the domain 0 and 
let O q be the unit square [0,1] x [0,1] . Assume that €H^(Q o > , 

i = l,...,n are given and denote by Z the linear span of , i = l,...n. 

Definition 7.1 Let p > 0 , e > 0 . By ^(Z.p.e.Q) we denote the 
family of functions [,€ H°(Q) such that for any square S = [a^,a^+h] x 
[a 2 »a 2 +h]CQ the following properties are fulfilled 

a) C€H 1 (S) 


6) There exists = z(h ^(x-a)), z €1 and constant M (both depending 


on C ) such that 


i) | U-Cj | 5 M-h 

H (S) 


ii) k-U , £ m 

H L (S) 


iii) k | 1 > pMh 

H (S) 


Let us illustrate our definition by a few examples. 

(i) Let Z be the set of all polynomials of degree less than or equal 

to n . Then any polynomial of degree £ n on Q belongs to the family 
H(Z,p,k,Q) with p and e arbitrary. 

ii) Let Z be the set of guadratic polynomials. Then £ = sin x^ 

belongs to the family H(Z,p,e,Q) for e £ 1 and some suitably chosen P . 
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iii) The family H is typically characterized by Z being the poly¬ 
nomials of degree <_ m , say, and then we take to be a suitable Taylor 

expansion of f . Then (i) and (ii) are more or less standard, and (iii) 
states that is not 'degenerate". 


LEMMA 7. 1. Suppose that f €t (Z,p,Q) then there exists C (depen¬ 
den t on (Z,p,e,Q), such that for any square SCQ 


(7.1) |f| < C(diam S)' 

n (S) 


inf 

d = constant 
functions on S 


1 0 
H (S) 


Proof . 1) Define F€H (Q q ) by 


F(x) = f(a + hx) 


Now by the assumption there exists F (x)€Z such that 

o 


(7.2) |IF-F I[ < M , 

H°(Q o ) 


(7.3) F-F , < M 

°«<v 


(7.4) F, > pMh 

° H 1 (Q ) “ 
o 


Denote now by F respectively F^ the average of F respective to F q on 

Q . Then ve have 
o 


(7.5) | | (F-F )-(F-F ) | i lIlF-Fjl <M 

° ° H°(Q ) ° H°(Q ) “ 

o o 
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because F-F 0 is t,le H ° projection onto the set of constant functions on 

% • 

Because Z is finite dimensional space there exists C (dependent on 
Z ) such that 


(76) If I = If -f I , < c ||f -f || 

0 H 1 (Q 0 ) ° ° H 1 (Q o ) ° ° H°(Q o ) 


and hence combining (7.4), (7.5) and (7.6) we get 


M < 


,£ -1 | 
h p 


° H 1 (Q q ) 


< Ch 


: P _ 1 |l F „-F 


H°(Q o ) 


and 


l|0 r -F o )-( F -F o ] 


H°(Q o ) 


Cp -1 h e 


I If -f ] 

J J o o * 


h ° ( Q q ) 


On the other hand 


F_F H°(Q ) " F °" F ° H°(Q ) " H U (QJ 

o o 


I|(F-F o )-( F - F 0 ) || 


and hence for Cp ^h £ < 1/2 we have 


(7.7) I IF-F|| 


> ^ F -F 

H°(Q 0 ) “ 2 ° ° n°(q 0 ) 


Further from (7.3) and (7.4) 
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(7.8) [F| < |F | + |F -F| . < |F | . 

H X (Q ) ° H 1 (Q ) ° H 1 (Q ) ° H 2 (Q ) 

o o o o 


+ M < IF i (1 + p’V) < 21F I 

H i (Q ) 0 H 1 (Q ) 

o o 


for h ( 'p 1 


Therefore by (7.6), (7.7) and (7.8), we have 


i < 2 I F I , < 2C||F -F I 

,1/^ X - 1 o' - 11 o O 1 


I| < 4C| |F-F| | 

H"(Q 0 ) " H j (Q o ) " ” ° h°(Q 0 ) ~ H°(Q o ) 


so upon rescaling back to S (7.1) is proven in the case diam S £ h o (p,e,Z) 
2) Suppose now that diam S > h (p,e,Z). Put 


r ,. _-iINT 
[diam S , . 

0 ' hr~J +1 


1 NT 

where [*] denotes as before the integral part. Clearly we can divide 

2 2 

S into a congruent squares S^, i = l,...,o with 


,. c. diam S , 

diam S. = —-—— < h 
i a — o 


Thus 


|£|\ - l 

H (S) 1=1 


2 

,|2 <c(^fLS) l inf | |f-d | | 2 

H (S ) i=l H°(S.) 

o ,_. . i 


d=constant 
function 


< Co 2 (diam S) 2 


inf ||f-d| 
d=constant 
function 


H°(S) 
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But o < 1 + —-^ < C and the result follows. 

— n ^ 

o 

k k 

We introduce the family ^(Z,P,e,Q) for Q = with 0~ as in section 

b b 

k 

2. Because fi consists of a finite number of Q. it is clear that we can 

b 

k 

extend the family H into H(Z,p,e,fl) so that the restriction on Q = Q_ 

b 

is the family fl(Z,p,e,fi) . 

In what follows we will assume that we are concerned with problems P(H o> w,g) 

introduced in section 6 where g H (g^,g 2 ) , g^ €H (Z ,p ,e ,0) . 

We shall discuss now the error estimate of section 6 in more detail. We 

have shown in theorem 6.1 that the essential part of the estimate is the norm 

of n which is defined on the star u> . In section 4 we introduced the 
P P 

standard star 4> (see lemma 4.3) as the image of to under the mapping 

J . 

P 

On $ we will define now the space H = H, * H„ with 
P P l.P 2, p 

H (* ) - tv - (v, ,v ? ) €H 1 (4> n ) x H 1 (4> n )|voJ €H(w )} 

PP L £ P ppp 


^tb ft(u'p) defined in (6.7). 

In 6.1 we defined the bilinear form B 


Let us define the form B 


defined on x with the same expression as in (6.1) but with integra- 

* - _ ^ ^ 
tion over , and A, p instead of A, y, A = A<>J , m = U°J 


LEMMA 7.2. The bilinear form B is such that for any u€H , 
---- ■ ___ __P___ P 


(7.9) C | |u | | 2 < B (u,u) < C 2 | |u | | 2 

1 n (♦ ) P HM* ) 


with constants , i = 1,2 dependent only on A,p and K but independent 
of P . 
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Proof. For each P€R(D) there are line segments T'’ C. 9<J>^ (i = 1,2) 

both being edges of standard P elements of <5>^ and such that u € Up = 0 

on T (i) . Applying now lemma 4.3 we see that there are only finite number 

of different cases of the domains <3> and the line segments 1'^'*'^ for each 

P 

the Korn inequality holds. Therefore (7.9) holds with | - | . replacing 

»<V 

| |■|| . . Because these two norms are equivalent with constants depending 

h (y 

on (y r (1) .F^ 2 -*) we get (7.9) immediately when using again lemma 4.3. 

Denote now by M (<t> ) the set of all u€H ($ ) being bilinear on 
P P P P 

every standard P element of the standard star <J> . Further we denote 

P 


e = e°J 1 , u = uojf 1 , U = U»J 1 , gj. = g^-T 1 , g £ = ^(diam A) 2 


i = 1,2 


where e, u, U, g were defined in section 6. 


We have 


B (e v) = 0 , Vv (<J> ) 
P P P 


8 p ( e ,v) = B p (U,v) - B p (u,v) 


*(“'V) = l )» Vv€H F 
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By integration by parts we get 


Bjl’.v) = l l 


[l 


A£5 i=l j = l 
P 


a. ,(U)n.v.ds 

1 , j 1 j 


. i (U)v.dxj 


3A A 

where (n^,n 2 ) is the unit normal on each edge of 3A pointing outward from 
A and 


„ 9U 3U 
a i,j = X( 9^ + 

„ , „ , 3 2 U- 


f 30 3U. I 

S 7 + ^J 

u I 1 J 


Li<U) - (A+m) 

. . . . sV 

4 <u > ■ n^sr 2 


Let F (respective r 9 ) be the union of all vertical ^horizontal) 

1 , P ^ » P 


edges of all A €. , Then 


(7.10) B (e,v) = l f-(fg.+L (U)],v ) 
P i=l L 


W + j-i (} ‘ <0 ‘> > -h ) .h<L, P >] 


where indicates the jump in a function across . The jump has 

an obvious sense if the relevant edge T of A is inside 4> . If 

P 

rC 3$ we have to distinguish the case of whether v.€H, =)v, = 0 on T 

P J 1»P j 

or not. If all v^ = 0 on T we will take the jump to be zero. If it is 

not the case, then is set equal to zero outside ♦ , and the jump 

is taken the usual sense. This convention will be used through the paper. 

Because U is bilinear on every A 6$ and are constant on A 

P 


we see 


that L^(U) are constant on every A and J^(o^ j(U)) are linear on 
every edge of A . 




Putting i’ ~ n ®J (6.8) gives 
P P P 


(7.11) B (ii ,v) = B (e,v) V v£H (<J> ) 
P P P P P 


Using Lemma 7.2 and proceeding in the standard fashion it is easily seen 
that n is the unique H (<t> ) function satisfying (7.11). 

p p p } ° 


LEMMA 7.3. There exists a constant C dependent only on K and p , A 
such that 


(7.12) 1|n 


P 1 
P H (t ) 


2 

1 C f l (||g 1 +L 1 (U))|| 

L i=l H 


(*p) 


2 

+ l 

j=l 


11 ^( 0 ^( 11)11 


Ar.)’] 


P roof . The set of possible S^'s is finite. Therefore we have by Sobolev 
imbedding theorem 


(7.13) ||v|I I c lIv|| 

H (r.) H l (4> ) 

P 

with C depending only on K . (7.13) with the Lemma 7.2 and Schwarz's Inequality 

leads to the desired result. 


LEMMA 7.4 . Assume that € H (Z, p, c , ft) then 


(7.14) 


| > C £ (ilgi+L (U)|| + l ||j.(d (U))|| ) 

H ■> ) i=l H°(4> ) j = l 1 1J H (T i ) 


where C depends only on K , the family U and A, p 








Proof. 

Suppose that the 

constant C does not exist. 

each n 

= 1,2,3, ... 



i) 

a K-mesh 



ii) 

p[n ] €R( p[n]) 



iii) 

i(nie " 1<4 p ( „]> 

with U 

^ being bilinear on 

iv) 


with 

!ll n] ! ! 1 c, ||[ n] - 

1 H (A) 


any A €5 with C 

pin] 

average of g[ n ^ on A 

independent of n and 


C Pin] 

I f< 

H°(A) 


such that the unique solution €H r ..(4 1 r ,) of 

pin] Pin] 

2 

r i *- . r i r 1 


(7.15) B 




+ j <J i <0 i,j <U » ,V j»-., In] ' ^Vnl'W • 

J-l L 2 ^ i 


satisfies 


(7.16) |h lnl || <± l (Iii| nl + L (U rn] )|| 


H it , ,) " 1-1 

P 1 " 1 


“ <* tn]’ 

P 


+ l IU,(»1 > I " , »II> „ 

j=i 1 ,J H°(r j ) 

l 


Putting r^ nl = |^ n] + L i (U Cnl ) 
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£.!”! = J. (O. , (U ln! )) 

1.1 1 l.J 


then without loss of generality we can assume 


(7.17) Y 
i = l 


ii „ * l iu! n !n „ - 


r ‘ 


h (r. 

i 


and 


(7.18) |r [ . n] ! . < c| |r* - M ni | | 

1 V(A) - 1 1 H°(A) 


VA €S 


n] 


The inequality (7.18) follows from the assumption that €H(Z,p ,e,Q) and that 
L^(U^ n ^) is constant on A . From (7.17) we have 


(7.19) i|r' 


A- r -[n]|| 

1 H°(A) 


1 I i r | 


H° (A) 


1 11 r i 


< l 


H <v 


Using Lemma 4.3 we may assume that <J> 


p(n) 


* ’ S [n] = S > H ^(0 fm1 ) = H 


J n] Jn]' 

P P 


Further by Rellich's lemma,(7.18) and (7.19), we may also assume that 


r. .t. in H°( I) . 

t i 


Now 


i . j i s 


i.J 


yieIds 


is linear on each edge of A€S and so we may assume also that 
H°(!') where T is the union of all edges of ACS. (7.17) 


(7.20) j ( ! |r 


i = l 


I o + IW^i ,H o )= 1 
H°(']’) i = l 1,J H (f) 


(7.15), (7.16) and (7.20) leads now to tlie contradiction. 
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LEMMA 7.5 . Let p€R(P) be such that = (v^, v^) € H(a)p) (see(6.7)) 


and u CQ. for some k€Z . Then 
_ p _o_ 


i In 


P u 1 


2 
C l 


H (<S> ) 1=1 

P 


+ I MViiWlI 


I ®i * * o L 1 1 ~ i~i \'' ' ' 1 o 

h (* ) j=i 1 1,J h (r.) 


A ; 


where restricted to any AC$ is the average value of g^ on A . 


Proof. The proof will be by contradiction. Suppose the lemma does not 
hold. Then for each n = 1,2,3,... we can find P^ n ^, pf n ^ , and 

gj n ^ as in i), iii), and iv) of the proof of lemma 7.4 when ii) is replaced 
by ii') 


ii') P^ €R(P n ) , (v ^,v ^ j)€H(w [ n ])»w^ n ] C Qe for some k€Z 


and (because of (7.10) and 15^(4,v) = 0) 


(7.21) 


-<i tnl + £ ± <0t")). v .j-J > o 

1 p ln] p lnJ H°(* [n] ) 

P l 


,[n], N _ T -1 


+ j I i < Ji < 0 iij( u"", ) ,v Jn] .rj n) ) u0 


= 0 


p' J p l J h (r ) 

i J n 1 
f»p 


[nj 


and with as in (7.15) and 


(7.22) | I n 


[n] 


> n 


£[n] A2[n) 


H ($ P [n]> Ll=1 


I (||g} nJ - 8 


H (* r _ 1 ) 


P t n J 


+ 1 IIV°i j (’J [nl ))ri 0 

j=i iJ M H°(r )J 


i.P 


[n]' 
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As in the proof of the previous theorem, by use of lemma 4.3 we can assume 
that 


a) S r i = S , <J> r , = 4> , T. , ,= T. , H r , = H are independent 

[n] pfn] i,pfn] l pin] 

P 


of n 


B) 


(1 , r C Q for some k£z where k is independent of n (because 


‘pfnl <a 
Z Q is finite) 


Now a), 6) yield that also Bp[ n j('»’) = B(*,*) does not depend on n 


') h 


fn) 


i i 

H ($) 


= 1 


Let us set now 




r< n h »' nl + LjCU 1 " 1 ) 


, - j.(o. ,(u (nl )) 

1,J J- 1,1 


Now by (7.21) and y we have 


| r [n] _ A—[nI 


1 H°(A) 


0 , VA€ $ 


(7.23) 


iu. .it - o 

X ’ J h (r.) 


I \ 

f 


\ 

l — 


It is easily demonstrated that 


(7.24) l ||L (U [n] ) - ? L (U [n] ) M < l IU || 

i=l H°(4>) i , j*=l H (T.) 



and so we must have 


(7.25) j |l.. (l' |n| ) - V L. (V ln| ) || ■- 0 

V(t) 


By (7.22), using lemma 7.1 Rellich's lemma and the fact that S is finite, 
we mav assume that 


(7.26) r[ nl r. in H°(l) 


(7.27) gj" 1 - g. Ln 1!°(;) 


Because of lemma 7.1 wo see that g. 


also in H^^(l') where 


H ($) is the usual space of fractional derivatives. Now by (7.23) again we 
see that r^ and g ; must be constant on each A€S . Further by lemma 4.3(i) 
and the fact that € H^* (1>) we see that g^. is in fact constant over all of 


By (7.26) and (7.27), it follows that 1.^(1’^) converges In H°(<t>) and by 
(7.25) it must converge to a constant on ■> . This follows that r^ is a constant 


By (7.21) in the limit we see that 


/v .-1 

(r.,v . «J . ) =0 

1 Ptn] Pfn 1 


Bv lemma 3.5 we have v . ,> 0 and v , 0 on a set of positive measure. 

pin]- pM 

Thus we conclude that r^ = 0 . 
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Returning to 7.15, we see that 


in! 


in H 1 ( 4-) 


which contradicts Y • 


Now we return to the question of" the a-posteriori estimate. fiiven a 
mesh V and the corresponding finite element solution l 1 ns in (6.5), 
we will associate to every AC'l the error indicator n(A) defined by 


(7.28) n 2 (A) 


l (diam A f | |J o (U) ||“ + (diam A) 2 ||g +L (U)|| 2 } 

i=i j=i 1 1J n O.a) 1 1 h° (a) 


where 


L^U) = (l+„) 
L 2 (U) = (A+u) 


o . . 

i, J 


3U 


1 


A ( 7hT + 

x i 



3 Xl 3x 2 


3x^3x 2 


3U. 

3U . 3U. 

uk—- + 

3x. 3x. 

1 i 

T 1 ) + 
3x 2 


and (J 2 ) indicates the jump across the vertical edges, 3^A (horizontal 

edges, 3 2 A ) of A . We use the same convention as before, i.e. if the edge 
rC3fi i s such that v.€H. -*■ v. = 0 on T then the lump J. is not taken 

ill J 1 i 

into consideration. If it is not the case, we take the iump J.(d. .) = a. . 

i i.l 

The estimator E then is defined 


E 2 = l r,“(A) 

A €0 

THEOREM 7.6 . T he estimator E is an upper and lower e stim a tor, i.e. there 
exists c^, c 9 dependent only on K, \, ij. Si and f but not u and V , such 
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that the error of the finite element solution of the problem P(H Q ,w,g) 
w€[M(P)]^ , g (Z t p ,t ,fl) satisfies the estimate 

(7.30) c.E < ||e|| 1 c, E 

1 1 _ 


Before proving the theorem let us prove two simple lemma. 

LEMMA 7.7 . Suppose F is an edge of some A£P . Then 

a) Either (i) FC3ft , 

or (ii) 3 Pf R(P) 

such that ACto but F(i3u) 

P ~ P 

g) r fl A ’ ^ 0 for at m ost x elements A 1 617 where x depends only 
on K . 


P roof. 

a) By lemma 3.5 we have Ev^ = 1 , so certainly for some P €R(P) 

v > 0 at the midpoint of T . Since v is a non negative and linear on T, 

P P 

then Vp > 0 on f except perhaps at an endpoint. The result is then immediate 
upon noting that v is continuous and inducing lemma 3.6. 

p 

0) The number of such elements is clearly bounded by the number with non 
empty intersection with A . Lemma 4.1 iii) then gives the result 

LEMMA 7.8. Suppose aj,...,a m are non negative real numbers and let a 
be a mapping from {l,...,n} -» (l,...,m) such that 
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i) a is ont o 

ii) 3 M ' s 0 such that Vl i j 1 m the set {i 11 ^ i _< n, a(i) = j} has 
at most M elements then 


m 

l 

j = l 


a. 

J 


n 

(i) 


m 

1 M l 
j = l 


a . 
J 


The lemma is obvious. 


Proof of the Theorem 7.6. Using a scaling argument applied to lemma 7.3 and 
7.4 we obtain 


(7.31) 



ii <... p > 


L 1 [(diam) 2 ||g +L (U) 
A€D i = l 11 

ACio 

P 


H°(A) 


+ diam A l ||J?(o (U ))|| 2 ] 

.1 = 1 1 1J H°(d.A) 


where J° indicates that we use the jump across edges 
but within L! equal to zero. Adding over all pCR(P) 
and lemma 4.1, we are in a situation covered by lemma 7 
Before formulating the next theorem, we prove 


A which are on 3w 


P 

and using lemma 7.8 
7 and the result follows. 


LEMMA 7.10. Let' R*(V) 


{P €R(i?) j v 

P 


0 on U 3qh, Then 

k€' Z 9 

O 


(7.32) U ‘i. ® Q 

P€R*(V) p 
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Proof. The lemma follows from lemma 3.2 applied to each separately 

-- Q 

which shows that any A€V has at least one vertex P say, being a regular node 

k k 

located in the interior of a . Clearly the corresponding , and so 

P€R*(V) . 

Let us introduce now another error indicator n(A) , 


(7.33) n 2 (A) = £ {diam A I||J.a (U) | | 2 + (diam A) 2 1 | g.- A gJ 1> 

n -*■ * J .-.Ox- , \ J- i ..O / . v 


«.= ! j=l 1 1,J H°O.A) 


H (A) 


where g^ is the average value of g^ on and the corresponding error 


estimator 


(7.34) E = l n (A) 

&€D 1 


THEOREM 7.11. The estimator (7.34) is also an upper and lower estimator. 


Proof. Let p€R*(V ) . Then using lemma 7.5 and 7.4, we see that for any 


(7.35) 


I (diam A) 2 ||g +L (U)|| 2 < f ((diam A) 2 ||g - A g || 

i=l H°(A) A6 w i=l 1 1 H°(A) 

P 


+ diam A £ ] |J°o . (U) | | 2 ) 

j = l 1 1,J H°O.A) 


We have shown in lemma 4.1 that any A is contained in not more than M(K) 
stars (i'p . Lemma 7.10 shows that every A is contained at least in one co^ 
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P fc R*(D) .md therefore summing (7.33) over all f\€V we see immediately that 

m — a. 

K < CK which yields that E is an upper estimator. Because E < E E is clearly 
a lower estimator. 


Lemma 7.12. There exist CL > 0 , i = 1,2 such that 


(7.36) C |jJ a . (U) 1 | 

H'().A) 


(A)m 2 


-t 

< diam A ) fJ.o • (U) (x^'') ] < 

ft=l 


c 2hVi,j <l 'dl'. 


HO.A) 

l 


(r\) 

where x ? are tlie vertices of A 


J.o 

l 


Proof. 


. . (!') is 

i, J 


The inequality (7.36) follows immediately from the fact that 
linear on every edge of A . 


Lemma 7.12 allows us to introduce another error indicator, 

2 2 4 

(7.37) n 2 (7) = ) {(diam A)“ ) j [J.o .(U)(x^)] 2 + (diam A) 2 || g . - A g .|| 2 } 

•■■■ l L j = l f.= l 1 1,1 * 1 1 H°(A) J 

and error estimator 

(7.38) E 2 = V p 2 (A) 

A6- ' 

and we have from lemma 7.12 and Theorem 7.11, 


THEOREM 7.13 . The estimator (7.38) is an upper and lower estimator. 
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8 . THE ASYMPTOTIC ESTIMATE 


We have shown in the previous section that the estimator E introduced in 
(7.35) and (7.34) is under the proper assumptions simultaneously an upper and 
lower estimator. 

In this section we will analyze an'estimator E , which will be equivalent 
with E i.e. there exist and so that 


(8.1) C 2 E 1 E < C t E , 


and will be asymptotically exact for the energy norm |||’||| introduced in (6.3). 
We shall say that an estimator E is asymptotically exact with respect to the 
energy norm if 


( 8 . 2 ) 


i -* 1 as 


To show (8.2) we have to make various assumptions about the solution u and 
the meshes in addition to the assumption that the mesh is a K-mesh. 

Suppose V(Q) is a mesh and let P'CP satisfy 

i) If A€P' , diam A = h , h > 0 

ii) if &€V then all vertices of A are proper nodes. 

If (i), and (ii) hold then we shall say that V' is uniform. If only 
(1) holds then we will say that V has uniform size. 

Suppose now that A C B C n then we shall write A < B if 
P(A,B) = dist(A,lR 2 -B) > 0 . 

Let P Q Cp(n) then we will write 
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( U ) - I n t 1 
■€D 


.ind ( L) ) M iiiibi i ■ >; el i-nients contained in V . We will write V" < V ' if 

o o 


,(p") <r s 


Finally denote bv M(P;’ ) the set of functions of (M(P)] restricted 

O 

to and M(P,) C M(p, ’ ) the set of functions which are zero at . 

We shall Make use of the following version of Theorem 5.2 of 111] . 


THEOREM 8._]_. Let P' C P(t :) he a unif orm mesh. Assume that the bilinear 

form B defi ned in_ (6.1) h as cons tan t coefficients on fl(P') and let V" < V . 

Then if \ = ('. ^, \ 0 ) € I (' (P 1 ) ] " satisfies 


(8.3) 


.3 ( ; ,:)! A / 


i=l 1 ' H 1 (.,(£?'> > 


for all i € f.t(V, )) then 


U,-x| 


(8. A) K.! , li c A f inf i.' •< ! | + ,, (off)"'. QOO ’) V o 

1 h 1 cap") i = i|v € '((p .(P)) l,( ^ p,) { { )} H°(«(P’)) 


( : i(V' r )(O' )) 


Il c (: (P')) 


with C dependent only on the bilinear form (6.1) (and not on fi(P") , ft(V ’)) . 

Suppose now that P'CP has uniform mesh size and let P" < V' . Then 
we can del in. a difference operator I. which maps H^(U(P')) into H (fl(P")) 
so t lia t 


(T.u)(x> = 3h | ii ( x+he . ) -u (x-lie •) | 

t l 1 






Then if V < V it follows that for D. = -r— , 

1 ax. we have 
1 
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Proof. For any w€f.!(P,;'(l}")) wo have 


(8.7; B(l : - u,w) = B(u-uu,w) = ) 15, (n-iiu,w) 

where ilu€'l(P) is such that ( !u) (x) = u(x) for any x€K(P), and B^, is the 
restriction of B to A . 

On each “\€.V" , B.(u-iiu,w) is the sum of terms of the form 

'W . 

——■- (u .-Hu .) v — dx 
x 2 1 J " N i 

for = 1,2 ;;nd some constant o . Consider first the terms 

tw. 




r f is oisv r«> see that for u a quadratic function the above expression is 
equal, to zero. Therefore by standard argument we get 


( 8 . 8 ) 




(u 'n ) 


Ch 


'»! 3 ■ w i I i 

! 11 5 (A) H (A) 


with 0 


hot 


independent °f ’A and u,w . 
us considet now the term 


(8.9) 



(u . -1'u . ) 
I 1 


aw 


ax 


dx. 


k * 


2 . 


Obviously we need only analyse the case k = 1 









I'sing the notation as shown in Figure 8.1 we integrate by parts in (8.9) 



Figure 8.1. The notation of an element. 



( 8 . 10 ) 
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<1 w. 

Ren 1 i i ng t ii.it (u — i: lj .) is j continuous fund ion on 1 (D") and •-— Is 

.1 J 3x 

iw. 1 

continuous on all horizontal edges of V" a ml = 0 on the horizontal edges 

' X I 

of <D") wo got by adding (8.10) over all € V" 


( 8 . 11 ) 


:i i , ? 

- (u.-liu.) - - — dx| •" Cl/' r o. ; 

j ox, j J )x. ' - ' ) 1 ,,3 


. 0 ") 




il H 1 (.Q TV")) 


for k i i 


(8.8), (8.11) and (8.7) yield 


!', f I •- 


< 1, ( 


) ( 


; i i tp"n 


i 


) 


H 1 (o(D")) 


Applying now theor• • fc. I ) wo got i rom (8. If.) 


(8.: i) .i:.ti.; • t 

1 1 II ( CD"’ ) j 




l + (o(D") nCD'" ) i - ' 1 ' 11 i I o 

H (u(D")) ( L K J H (n(V" 


< . (£>’" ) . (O' 


•y ). { I,z i '"j i ! ) + I ! ! y _u ; ! i 0 j 

; = I 1 II (2?")) 1 1 H (P(0"))J 


whore wo have used ft,, t ad t ii.it 


(8.14) M -’a. ! 1 • j jl'.-u. | j + j ] u ,-llti | | 

1 1 ll‘'( ■XV") ) 1 1 H°(t.(P") ) 3 J H°(;?(p")) 


H°(,l(P’ 


Cl/ 


)) 


l"l I 


ir('.i(P"» 


an* 
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Let us consider further D.(uj-riu.) on A £0'" . Refering to Figure 8.2 



Figure 8.2. The notation of an element. 


2 A 

for any w£H (A) we choose 4>^ = so that 


(Z (l) ) = w(Z (l) ) 


where Z^^ are the midpoints of the particular sides as shown in Figure 8.2. 
For u. = g a quardatic function , we can easily check that 


And 


D^g-Hg) = J i D i g . 


so by the usual arguments we have that for some ,€H.(A) , 

3 • 1 1 


(8.15) ||D (u -flu )-♦* || <Ch 2 |u j 

1 3 J J>1 H (A) J H J (A) 


Summing over all &€V" 


we obtain 












(8.16) 


i). (u ,-n». ) 
i j i 


i > * ‘ 


Ch“]u.j 


u i .Kir' )) j n 3 (n(P"';) 


Combining now (8.16) and (8.13) wo ^et the desired result. 


hKMMA_ 6. A. L’sinj; i lie sane notation ;is i n T heorem 8.3 we get for k 


C’.Cu.-h. ) . h, in.-V. ) ) i - C 

1 1 w‘\ a v ")) 


■a 


11 ; I i a 

J H (i: (V")) 




i (r > i 




h" u 


!„ 1 ; 


+ ! u ; -U : 


■' '., 3 ( ‘(V")) J j 


wnere 


( ( L ! ’" , ! V") i 


1' roe i . Wr i 1 1 


( (O'” ) , (P")) ( V ii 


i=i 1 ii S (V")) 


4 U . 


1 ! 


I 11° f r 


ir (.<;(£>")) 


Now we have ior any € 0" 


(8.17) I'D i a (n.-il. ) ) 

• 1 1 k ’ 1 n n (M 


= (;•: . +i>, (u .-r.) - 

I . ' 5 3 1 


'':.'.'^i,k + °k (l V 


, , » - < . ) + (n,(n.-u.) - <^ A „,<j> a > 

’’ K H°C,) ' H "(A) 1 1 » J ’ ? * ,k ll°(A) 


4- (;.'. )-< «»i). (u.-u,)~ < f>, ,) 


3 LA 


1 

°(n(P"’))J 








= 0 . Summing (8.17) 


It is easy to see that for S' ^ k we have (<j>^ „»<(>. ,) 

i,h H o (A) 

over all a€V" and using Theorem 8.3 we get 


(8.18) | (D (u -U ),D (u -U ))! 

' J J k 1 1 H°( fi (^ 1 " i) 


<C((» 2 -H.( l l u* I 2 ) 1/2 ) 

i,j = l A€P"' ,J H°(A) 


Using once more Theorem 8.3, we get 


(8.19) 


v A 

l « 


J ’ 1 u° 


tT(fi(P ,M )) 


< C[<H-| |D (u -U ) | | 


i J J u° 


H u (n(P"’ )) 


Using now Theorem 8.i for £ = u-U we have 


(8.20) ||D (u-U .)!! 

1 J J H°(fl(P'”)) 


4 . 

£ C( l (h||u || + 

j = l J IT (£!(!?")) 


3 3 H°(n(P")) 


and combining now (8.18), (8.19), and (8.20), we obtain the desired result. 


LEMMA 8.5. Let S be the square as in Figure 8.1 with side length h 
Then for any 


f = ax^ + bx£ + cx.^2 + 8 


we have 


(8.21) | | ax. | j 2 - h 2 (f 2 (N + )+f 2 (N')+f 2 (S~)+f 2 (S + ))| 

1 H°(S) 


< M ||bx 9 +cx.x 9 +d j I 

H°(S) 


with M independent of h, a, b, c, d. 
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Proof. We have 


(8.24) D.(u,-U.) = (D.u.-T h U.) + (T h U.-D.U.) 

iJJ ijij ijij 


By Theorem 8.2 we get 


(8.25) 1 |D u - tJu || 1C l l (h 2 ||u || + 

J 1J H 0 (U(P"' )) k=l H J (fi(P')) 


+ [p(fi(P"' ),ft(P"))]~ 1 | |u -U | ! ] 

R R H°(n(P')) 


Theorem 8.3 yields 


(8.26) ||D( U -U)- l / || 

J j a€P"' J>1 h (a(P"’)) 


1 ciD _1 (n(P" , ),n(P")) I 


k=l 


{h2 !l\!l 3 + llu k -U || >1 = E 

H J (w(P") k H°(n(P")) 1 


Combining (8.24), (8.25) and (8.26) we get 


(8.27) 


j |t'V -D H.- 7 *■' || 

1 J i i » l u° 


< C E, 


, €J ", H w (fi(P'” )) 


Simple computation shows that T.U.-DjU^ is bilinear on each A€P"' . For sake 

of definitions and without loss of any generality let us suppose that i = 1. 

Consider now T^Uj-D^Uj . Using notation shown in Figure 8.3 we have 
for l = ( f , 1 ,^ 2 )€r E 


r 
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N' NT 



Figure 8.3. The scheme of the notation. 


(8.28) (Tju -D^HO 


2h 2 [U (F.j+h,?^ - U.(C r h,? 2 )] - h _1 [Uj (F 1 “h,C 2 ) ] 


1 

2 



For we get analogously 

(8.29) (tV-D U )(c) = - \ J. AO 

h 

Let now <J> be the L„ projection of (T U.-D U.) onto H (A) . Then from 

^ * J J -I* 

(8.27) we get 


(8.30) 


IVj-Vj) 


h (n(P"')) 


: I 

A€V"' 


iK^ 2 o +o(e i> 

H°(A) 


Now it is easy to see that for any bilinear function f on A , f * ax^ + 
bx 2 + cXjXj + d (with coordinate origin in the center of A ) the L 2 projec¬ 
tion of f onto H^(A) is exactly ax^ and so by lemma 8.5, (8.28) and (8.29) 
we obtain 






(8.31) l ||$ M = fg l fJ (N“) + J (N + ) + J (S ) + J 

/\£P"< U°i'A'V 48 • A^n»" i »3 l 


H (A) 


A€P 


0(||tV-DU- l $Jj| ) • 

J 33 A€V'" 3,3 H°(Q(P"' )) 


using now (8.27) and (8.30) we get 


(8.32) HtV-DU || 

3 3 1 H°( P(P"' )) 


<l>? ,(P"’ ) + 0(F. 2 ) 
A * J A 


Hence we have 


(8.35) | (D (u-U)|| 2 "'~ h 

1 J 3 H°(n(P’") 


| | (T*?U .-D,U.) + (D u .-T^U . ) ] | 

1 j 1 j i j i j 1 Vw") 


Irt'u.- d,u. ' |2 


13 ^ H<W)) 


+ HD u -T^U || 2 + 

13 L 3 H (.HP'" ) 


+ «||D u -T U || I|T"U -D U || 

13 33 H 0 (ij(P"')) 3 J H°(n(P'")) 


where -2 < n < 2 

Further using (8.25) and (8.20) 


(8.34) ||tV-du|| < I |d (u -U ) 11 + 

3 3 33 H°(n(p")) 3 3 3 H (n(P'")) 


+ T 


JVVill O - C(E 1 + E h fl u ill 2 

3 3 3 H°(n(p")) j=i 3 ir 


(n(P’)) 


+ p _ 1 (fi(P" 


| U . -II . | | 

3 3 H°(S?tP')) 


,(S + )) + 

9 J 


>,n(t?")) 


Using now (8.32), (8.33), (8.34) and (8.25), we have 
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(8.35) 


l D i (u r u i ) M 2 o = *1 + 0(E ? + E i ( l MM 3 + 

1 J 3 H°(D") 1>J 1 1 j=l J H 3 (n(t>')) 


+ p" 1 (fi(P" , ) f n(p"))||u.-u,|] ) 


3 3 H°(fi(P')) 


and (8.35) yields almost immediately the theorem. 

In section 6 and 7 we introduced various error indicators. See e.g. 

theorems 7.11 and 7.12. It is obvious that we can have many equivalent error 

indicators which would be simultaneously upper and lower ones. Theorem 8.6 

enables us to design a special one, optimally suited to our purpose. 

In (6.1) we introduced the basic bilinear form. Obviously we can write 
T 

B(u,v) = [Dv]A[Du] where 


[Du] - [D 1 u 1 ,D 1 u 2 ,D 2 u 1 ,D 2 u 2 ] 


and analogously Dv 

The matrix A has then the form. 


- 



- 


— 


X+2u 

0 

0 

X 




0 

u 

u 

0 


A 11 

A 12 

0 

u 

u 

0 




X 



X+2p 


A 21 

CM 

CM 

< 

_ 



_ 


_ 

J 


Assume now that A€V'" as in theorem 8.2. Then define 


Q = j£ 1 {fJ l,l’ J l,2 jA ll t - 7 l,l’ J l,2^ [a j ] + ^ J 2,1 ,J 2,2^ A 22 [J 2,1 ,J 2,2 ] 


[a J 1} 
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where a. j = 1,2,3,A 
Further let 


are the four vertices of the element 


A . 


B‘(A) 


1 


i=l 



2 


H° (A) 


and then define the error indicator (which obviously is simultaneously an upper 
and lower one) 


(8.36) n 2 (A) = |g a 2 (A) + h 2 - Y 6 2 (A) 


where y>0 is a constant which will be determined later. We mention here only 
that for smooth g we have J (3 ( A ) j < Ch . 

Now we have 


THEOREM 8.7. Let the assumptions of Theorem 8.6 be satisfied. Then 


(8.37) HIu-U| 


.1(17" ) 


= l a 2 (A) + C[ l {(p V+p’VHtu, " 2 


usy' 


2 

l 

k=i 


k 1 1 ^ 

H J (fi(0')) 


+ P'’ 2 |l e l! 2 + (p 1 h+p“ 2 h 2 ) | |u,| | 3 lie | | 

H° * H J (Q(P’)) H° 


+ { l h l|u.|| 3 

j=l J H J (n(0')) 


+ p 1 \ |e|| 1 < l h i |U || + 

H°(n(P’)) j=i J H J (n(P')) 


+ e 


} 


H u (ft(P’)) 


Proof. We have for u - U = e = (e^,e 2 ) 
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j \ (A+2n)((D ie] ) 2 +(D 9 e 2 ) 2 ) + u[(U 2k s 1 ) 2 +(D ] e 2 ) 2 ] + 2p[D 2 e 1 D 1 e 2 ] + 
fi(P’” ) 

+ 2\ [D^e^D 0 e 2 J}dx 


Using Lemma 8.4 and Theorem 8.6 we get the desired result. 


So far we have defined the error indicator n(A) by (8.36) for all A€ V” . 


We will extend this definition to all elements A with A€Q (see sections 

u 

2 and 6) where the coefficients are constant and A LI 3Q = 0 . If A has an 

U 

edge on 3Q then we shall define the indicator in terms of the jumps in o. . , 
” i»j 

so that we obtain an upper and lower estimator. The detailed formulation arid 
extensive numerical experience will be addressed in a forthcoming paper. 

Here we will assume only 

a) The indicator has the form (8.36) for any A 6, Afl3Q^ = 0 (y > 0) 
fixed but arbitrary. 

b) All the meshes on K-meshes and the indicator leads to a simultaneous 
upper and lower estimator (under the assumptions spelled out in Section 7 ). 

We will place some additional assumptions on the solution u and the meshes 
used, as listed below. Then we shall prove that the error estimator is asymp¬ 
totically correct. By this we mean that 



as || | e | |! -*• 0 . 

We will assume that there is a sequence of K-meshes 
h ■* 0 with K independent of h . 


h. 




:) The mesh is equilibrated in the sense that max n(A)/ min n(A) <_ C 


with C independent of h . 

d) !|)e|!| ;• Ch with C > 0 independent of h . 

e) There exists s > 0 and C , both independent of h , such that 


MIL ±a,|||ei 

h (n) 


f) For each e > 0 , li > 0 , there exists meshes V’ . > V" > V'” 

h,t,ih,e,ih,£, 

i = l,2...m(h,e) such that 

i) the bilinear form has constant coefficients on ST(l?' .) 

h,c ,i 

ii) V " . are uniform and V' . are of uniform size 

h,£ ,i h ,£ ,i 

iii) n(£>"' .)nfl(K" .) = 0 for all i ¥ j . 

n , l , l n, c , 1 

iv) p(ST(D" .), .)) > Ch 6 with C and 9 independent of 

h, t , l h, e , i — 

h and e , 0 < 0 < S 

m(h,£) 

g) Denote R) = U S2(P' .) and analogously R" and Rl*' and 

n, e i — 1 n ,£ ji n, e rijc. 

we shall assume that 

i) | I u | | 1 X(e) 

H J (R' ) 
h ,£- 

ii) | |g| | , i X(e) 

>• (R h. c > 


m(r,h) 6(h) 

h) Let S’" = 9. - U •) and S"' be the 6 neighborhoou 

h,£ h,e,i h , e 

e= l 

of S"' with * = h° 0 < 9 •- S . 

h,E 

Denote now 


& . = (A IA n S"’ ,Uh) t 0} , 
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flncf we will assume that 


o(i) ) <- ro (P) 

h,t - h 


where o(P^) respective o(ft^ f ) denotes the number of elements contained in P^ 


respective ft, 

h,f. 

Now we have 


THEOREM 8.8. Suppose that the assumptions listed above are satisfied, then 


e >0 such that for e < e and h < h (c) 
o__o_— o 


< Mr 


where M is independent of h and t: 


Proof. For every Pj|' ^ ^ we can write by theorem 8.7 and property f) and 


g) 


m(h,r) 0 

? ^ e ^fi(P"’ .) 

i=l h,E,i 


l n (A) + z - 2 


' € U V"h E i 

i = l 


where 


3-O..4-20,„2,., , L -20 + 2s |i|e(i| 2 + ^1+3-6^2+3-20, x(e) 0 |,| q + 


|z|< C[(h J ”+h™)x"(e) + h 


+ (h 1 +S - 6 X(r) + h ZS - Zt, |||e||| 0 )( X (E)h^ S +|||e|||)) 


2s-26 


2-s, 
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and 


| z| = 


■yh 


m 

A€U V"' 
i=l 


J I! g i' g i 

1=1 


h , r , i 


2 

H° (A) 


Taking into consideration the property d) we get 


z < C 


(h 1 - 0 4h 2 - 2s )x 2 (e) + h 2(s - 9) + 


(h s-e +h i + s-20 


/L l-0 2, v ^ . l+s-20 , N 
(h X (c) + h x(e) 


, uS-6 , x . . 2s-20, 
+ h x(e) + h ] 


We have 



2 

H°(A) 


< 


Ch 2 ||g.|| 


2 

H L (A) 


and therefore 


ch 1 I) 

i=l h ,c,i 


Ch 


X(e) 


Hence 


\z\ + \z\± 


,2 

S2 


for h < h^(f) 


) x (e) + 


We have further 
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J. n 2 (A) 

A€P t 

h 


and by property 


2 

n 

max 

2 

n A 

min 


1 C cl n 2 (A) 
A€D 

h 


when using property c) . 

So we have for sufficiently small e 


l n 2 (A) < o(b )n 2 (A) < r.o(P )n 2 < e £ n 2 (A) 

£ h,t max - h max - in 

\€V. At h 


m 

\f u n '". 

r . ,1/ h,e,i 
i=l 


^ n (A) + , l n (A) 


m 

,\€ U V' 

i = 1 


h , r , i 


MED. 

n 

. m 

A f U P"' . 

4-1 h,t,i 


h) we get 


(1-eM) )' n 2 (A) 

A€P U 


n 2 (A) 


m 

A€ U P'" 

i=l 


h, c, i 


£ y n 2 (A) 
A€P 

h 


Using further the .act that the estimator is an upper and lower estimator we get 


(».39) y 11 Ie 1 I I$207*’* ) = E ‘ (1+P(h,e)r) 

f = 1 V h , f , i ’ 


(h, ) j ■ M , independent of h and e 

i I v we have 


m(h 

) 

\ - 


I 2 

'(KP’" 


h , r , i 




|2 

'S'" 


h ,e 




—— 


i 2 

Let us analyze now |||e|||g. 


It is easv to see that there exists a function 


^ such that 


i) 0 < iJj. < 1 on 

h — 

“> l D \! ;- s ThTi Ch ‘° 


iii) « h = 1 on s- 


Ip, = 0 on ft - s"' 6(h) 
h h,c 


Obviously now 


(8.40) III e ||| s ,„ £ B(eiJ) h ,eti» h ) 
h,e 


A typical term of B(e ^ is 


A D 1 (e.i|/ h )D k (e (l i(J h )dx 


where A is piecewise constant. 


We have now 


D 1 (e^ h )D k (e f il h ) = D 1 e^D i (tj^ep > + z* 


where 


z * = V^hVS + ^h Dl Vj D %. + e j e * D V\ - 2 v\° le j e £ 


Az*dx | £ c|j|e||| 2 [h s 0 +h 2s 20 ] 


Hence 
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and so wo have 


(8.41) B(ei^ h ,ev h ) = B(e,e'|>~) + R 


where 


(8.42) R < e|||u||| 


, 2 


if h £ h Q (e) 


Let us analyze now the term Btejep^) . Obviously supp et^CSKD^ £ ) . 

Using theorem 5.9 and the remark to it and lemma 4.1, we can find w€’M(P) , 

% 

r \j 

supp wCfi(P h £ ) such that 


2 

l I |v (e^-w) || 2 < c| | e'lv | I 

V€R(V) p h H 1 (ii) h H 1 (f2) 


O/ 


% 


and a(V ) < Ko(P ) 
h, e — h , f. 

Denote further 


«S.« ' <*•«*»,,> l-pOBW,, ) -«> . 


Then repeating arguments of theorem 6.1 we get 


(8.43) | B(e, ip. e) | <_ C 


2., J l 


r 1 iii" P i | i 2 i 

LP€R* p J 

h, : 


2| 1/2 2 
n p ll| | III <l' h e | 


Denote now 
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V. = {A | A C w , P€R* } 
h, g 1 p h, c 

Using lemma 4*1 we easily get 

_ •x, 

(8.44) a(V u ) < Ka (V ) < K a (D u ) 
h,e — h,e — o h,f 


and 


(8.45) l 

P ^, ( 


! 2 < C n 2 (A) 


A€P, 


h ,e 


In fact if on w the bilinear form has constant coefficients and w 03f2 = 0 
P p 

then lemma 7.5 shows that III Hi can be estimated from above by the suitable 
turn of the error indicators. In the other cases we first use lemma - 7.3 and 
lemma 7.5 and the same argument as used in theorem 7.8. Therefore we get 


|B(e,if; 2 e) | < Cf n 2 (A) ] 1/2 | | ifrje | | 


ACJ 


h,e 


Now 


\^\ 


| < c IU 2 e|| lC[||^e|| 2 + 

H 1 (n) h h 1 (S'" ) 

n y £ 


'^"^>-S h > + 

h, e h ,e 


+ ll ^"„ 2 w 5 <-v 

h,e 

because tj>, ~ 1 on RJ" and ip, = 0 on ft - S'" we get 

n n, e h h, e 


(8.46) ||U 2 e!|| fi < C(| IVe" 2 


,2 m2 


h e I I i + I!^ h e l I i 6fhl ) 

H i (S"' ) h H 1 (S'" (h) -S'" .) 

n j c n, g h j t 
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Consider 

(8.47) 

Now 

(8.48) 

and 

(8.49) 

(8.50) 

For x€S 


now 


'^H 1 (S'” ^ (h) -S)" ) 
n, n , r 


^ H ° {s h" e (h)-S h e ) ^ he ^H 1 (s;" ^ (h) -S" u J 


h,e h,e 


' ^ ^H°(S"’ 6(h) - S ’" ) 

h,e h,e 


^H 1 (S'" ^ (h) -S'" ) 
h,e h,c 


6(h) -s l " ) i=l 
h,e h,e 


= X[El I 


H A (S ,,f v - t ^ ^ 


H (S’" ~ V “ / -S"' ) 
h,e h,e 


+ | U,D 1 (^. e) || ... ] 

^ ^ U°/C»M (h) C ||» \ 


*T(S' M ) 

h, e h ,e 


l >, D X <0 e) | | 

n h »^ (h) pm • 


1 I |D i (i|), e) || . . 

H -(g,., V V „, .) h H°(S'" 6(h) -S"' ) 

ti|E n, e h,e h, e 


< I i^h e l I 1 C/.N 

h h 1 (S'" 5 (h) -S"' ) 
h, z h, e 


| I (D 1 ^, H. e| I 


a< h , <. c »- 9|| vII hU S (h). s „, 

h,c n, f. h,e h, e 


in ^ (h) o'n 
ti.e h,c 


we have 
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dist(x,S2-S'" * (h) ) < h G 
h, e 


and so 


IKell 


h 1 ' o 


H °^ S h' c (h)-S h' c ) ^(S"' 5(h) _ s ^'» ) 


h,e h.e 2 


o(h) 

J h,f. 

So we have from (8.43), (8.49) and (8.50), 


,2 ,2 


h ^(S"’ 6 (h) -s:" ) ^(S"' 6 (h) -S"' ) 

h,c h,e n h,e h,e' 


Recalling (8.46), (8.47), we get 




and therefore 


1/2 


(8.51) i C o(K h>e ) n ma J |* h e| 1^ < C E 1/2 (o«y ) V \ a J | V | | 


< Ce 


1/2 


f l n 2 <4) 

H. . 


1/2 


•V s 11 i, ! Ilvll i 

H (fi) H (fi) H (fi) 


where we used the fact that the mesh is equilibrated in the sense of the assumption 
c) above and that the estimator is an upper and lower one. 


Returning to (8.41) 








{Wit 
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ivii 2 i 

H A (ft) 


< B0|>. 


A e) 


Cr 1/2 ( 1 


i 6 '>H 1 («) + 


I I 

h (n) 


v 


» i * 


and so 


IVI l'i 

H A (n) 


< Ce 


1 1 


. e. 


B(ip. e,i|). e) 1 Ce | |e| | 2 1 

H (ft) 


and so by (8.40) 

IIWIlL I ce||e|| 2 . 

b h,e H 1 (ft) 

and the result follows. 

The theorem 8.8 obviously shows that the error estimator is asymptotically 
correct. This property is achieved for any y > 0 . Let us discuss now the 
selection of y . Obviously when the finite element solution is identically zero 
then the term in the error indicator associated with the jumps on the edges will 
disappear and the error indicator consists only of the "volume" integral 

, 2 o 

Y(diam A) [ ||g.-g || 

e=l 1 1 H°(A) 

2 

Assume now that we have a uniform mesh in IR and the exact solution is periodic 







u i = [sln IT x i Sln IT x 2 ]c i 

r . 2t\ 21 | 

u 2 = [sin x^ sin —- x 2 ]C 2 


Then it is easy to see that the finite element solution is identically zero. This 
leads to the choice 


(8.51) 


Y 


1_ X+3u 
!lT (X+2p) 2 +X 2 


As said above, the choice of y is rather arbitrary and (8.47) is one of the 
possibilities. 

We have assumed many very particular properties of the meshes and solutions, 
in the theorem 8.8. The problem arises whether these conditions can be satisfied. 
Practically we create the meshes in an adaptive mode. The experience has shown 
that the meshes which are adaptively constructed have roughly these properties, 
and that the effectivity index, seems to converge to 1 quicker than the 

theorem 8.8 supports. 

2 

For uniform meshes in IR and smooth solutions it can be shown that 


ll|e|f| Q - fc (l+0( E 2 )) 


(see [11]) and adaptively created meshes seem to lead to the same behaviour. 
In the next section we will discuss some concrete illustrative examples 


pertinent to these questions. 
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9. THE FEARS PROGRAM 


Based on the theory explained above, the program FEARS (Finite Element Adap¬ 
tive Research Solver) was developed. FEARS is a fully adaptive program solving 
a system of two elliptic equations and produces the error estimation (in 
various norms) together with the numerical solution of the given partial dif¬ 
ferential equation. The admissible domain is a union of curvilinear rectangles. 

The adaptive approach is based on equilibration of the error indicators. Ihe 
description and experimentation with FEARS will be reported elsewhere. 

In this paper we are using FEARS as an illustration of the developed theory. 

We will discuss here two examples. In both, we are concerned with the (plane 
strain) elasticity problem. We assume that E = 1 (E is the Young's elasticity 
modulus) and v = .3 (v is the Poission ratio). 

Example 1 . The elasticity problem on the square with displacements pre¬ 
scribed on the boundary. The data are shown in Figure 9.1. It is easy to see that 
the solution belongs to the space r (fi) (e > 0 , arbitrary). 



Figure 9.1. The data of example 1. 
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Example 2 . The elasticity problem on the rectangle with mixed boundary con¬ 
ditions. The data are shown in Figure 9.2. The solution is of the "stamp" type 
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Figure 9.2. The data of example 2. 

and the singularities are of the type described in [12] . Solution belong to 
3 / 2— e 

the space H (Q) , e > 0 arbitrary. 

In both cases the exact solution is not known, nevertheless by now elaborate 
computations we estimated with sufficient accuracy the exact energy of the solu¬ 
tion. This gives the possibility to compute the (exact) energy norm of the error 
and compare it with the estimator. 

Example 1 . Because of the obvious symmetries of the solution, we can compute 
the solution only on the quarter of the original square applying boundary condi¬ 


tions shown in Figure 9.3. 









Figure 9.3. The boundary conditions for example 1. 

FEARS constructs adaptively the meshes by equilibrating the error indicators. 

Figure 9 . Aa,b,c,d ,e , f, g show the sequence of constructed meshes. We see that 
the sequence of meshes satisfies the assumptions made in section 8. 


Figure 9.A. The sequence of adaptively constructed meshes for example 1. 
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Figure 9.5 shows the dependence of the energy norm of the error on the number 

of cne elements X . The norm at the error is measured in percent of the energy 

2 — € 

norm of the solution (|||u|||) . Because the solution belongs to H (ft) , 
the rate of convergence is N 2 (or more precisely N 2 ) for the uniform mesh 


as for the adaptive one. The rate N 2 is the maximal possible rate because of 
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We aisosee that the quality of the error estimator is better for the adaptively 
constructed meshes than for the uniform ones. This is likely the consequence 
of the equilibration of the error indicators which is essential in our theory. 

Example 2 . Because of the obvious symmetries of the solution, we can 
restrict the computation to the domain on a boundary condition shown in Figure 9.7. 










Figure 9.8a, b, c, d, e, f, g, h, i shows the sequence of the meshes constructed 
by FEARS. Once more we see that the assumption about the mesh made in section 8 
is essentially satisfied. 



Figure 9.8. The sequence of adaptively constructed meshes for example 2. 

Figure 9.9 shows the behavior of the energy norm of the error. Because the 

n / n "X / 0 

solution belongs to H ' (r > 0 arbitrary) and u H the rate of convergence 

_i. 

of the uniform mesh is N u . This is in complete agreement with the data shown 

-h 

in Figure 9.9. The adaptive mesh gives the rate of convergence N which is the 
maximal possible rate for the smooth solution. We see that the adaptive mesh 












































































Figure 9.9. The energy norm of the error — Example 2. 

removes the influence of the singularities on the rate of convergence. We see 
also very clearly that using a uniform mesh we practically can never achieve an 
accuracy of 5%. 

Figure 9.10 shows the behavior of the effectivity index for the Example 2. 


Once more we see that the effectivity index has practically acceptable value when 
the accuracy of the solution is in the range of 5-10%. In addition the rate of 
convergence of the effectivity index seems to be twice as high as that of the 
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Figure 9.10. The effectivity index — Example 2. 


solution. Also we see that the error estimator performs much better for the 
adaptively constructed meshes which equilibrate the error indicators than for 


the uniform mesh. 












103 


REFERENCES 

[1] P. Zave, G.E. Cole, "A Quantitative Evaluation of the Feasibility of, and 
Suitable Hardware Architectures for, an Adaptive, Parallel Finite-Element 
System". In preparation. 

[2] I. Babuska, W.C. Rheinboldt, "Reliable Error Estimation and Mesh Adaptation 
for the Finite Element Method", Computational Methods in Nonlinear Mechanics, 
J.T. Oden, ed., North Holland Pub. Co. (1980), pp.67-108. 

[3] P. Zave, W.C. Rheinboldt, "Design of an Adaptive Parallel Finite-Element 
System", ACM Trans, on Math. Software, Vol. 5 (1979), pp.1-17. 

[4] W.C. Rheinboldt, Ch. Mesztenyi, "On a Data Structure for Adaptive Finite 
Element Mesh Refinement", ACM Trans, on Math. Software, Vol. 6 (1980), 
pp.166-187. 

[5] I. Babuska, "Reliable A-posteriori Estimates and Adaptive Procedures for 
Solving Partial Differential Equations”, Finite Elements in Water Resources, 
S.Y. Wang, C.V. Alnso, G.F. Pinder, C.A. Brebbia, W.G. Gray, eds., School of 
Engineering, Univ. of Mississippi (1980). 

[6] I. Babuska, W.C. Rheinboldt, "A-posteriori Error Estimates for the Finite 
Element Method", Int. Journal Num. Mech. Eng., Vol. 12 (1978), pp.1397-1615. 

[7] I. Babuska, W.C. Rheinboldt, "Analyses of Optimal Finite-Element Meshes in 
R 1 ", Math, of Comp. Vol. 33 (1979), pp.435-463. 

[8] I. Babuska, W.C. Rheinboldt, "A-posteriori Error Analyses of Finite Element 
Solution for One-Dimensional Problems", SIAM J.Num.Anal. .18 (1981), pp.565-589. 

(91 I. Babuska, W.C. Rheinboldt, "Error Estimates for Adaptive Finite Element 
Computations", STAM J. Num. Anal. 15 (1978), pp.736-754. 

[10] J.A. Nitsche, A.H. Schatz, "Interior Estimation for Kitz-Galerkin Methods", 
Math, of Comp. Vol. 18 (1974), pp.937-958. 





104 


[11] A. Miller , Dissertation, (1981). 

[12] N.I. Muskhelishvili, "Some Basic Problems of the Mathematical Hieory of 
Elasticity", P. Noordhoff (1963), chapter 19. 






The Laboratory for Numerical Analysis is an integral part of the Institute 
for Physical Science and Technology of the University of Maryland, under the 
general administration of the Director, Institute for Physical Science and 
Technology. It has the following goals: 

• To conduct research in the mathematical theory and computational 
implementation of numerical analysis and related topics, with emphasis 
on the numerical treatment of linear and nonlinear differential equa¬ 
tions and problems in linear and nonlinear algebra. 

• To help bridge gaps between computational directions in engineering, 
physics, etc. and those in the mathematical community. 

• To provide a limited consulting service in all areas of numerical 
mathematics to the University as a whole, and also to government 
agencies and industries in the State of Maryland and the Washington 
Metropolitan area. 

• To assist with the education of numerical analysts, especially at the 
postdoctoral level, in conjunction with the Interdisciplinary Applied 
Mathematics Program and the programs of the Mathematics and Computer 
Science Departments. This includes active collaboration with government 
agencies such as the National Bureau of Standards. 

• To be an international center of study and research for foreign students 
in numerical mathematics who are supported by foreign governments or 
exchange agencies (Fulbright, etc.). 

Further information may be obtained from Professor I. Babu§ka, Chairman, 
Laboratory for Numerical Analysis, Institute for Physical Science and 
Technology, University of Maryland, College Park, Maryland 20742. 









ILME 


